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Abstract: Soft gluon exponentiation in non-abelian gauge theories can be described in 
terms of webs. So far this description has been restricted to amplitudes with two hard 
partons, where webs were defined as the colour-connected subset of diagrams. Here we 
generalise the concept of webs to the multi-leg case, where the hard interaction involves 
non-trivial colour flow. Using the replica trick from statistical physics we solve the combi- 
natorial problem of non-abelian exponentiation to all orders. In particular, we derive an 
algorithm for computing the colour factor associated with any given diagram in the expo- 
nent. The emerging result is exponentiation of a sum of webs, where each web is a linear 
combination of a subset of diagrams that are mutually related by permuting the eikonal 
gluon attachments to each hard parton. These linear combinations are responsible for 
partial cancellation of subdivergences, conforming with the renormalization of a multi-leg 
eikonal vertex. We also discuss the generalisation of exponentiation properties to beyond 
the eikonal approximation. 
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1. Introduction 



1.1 Wilson lines and exponentiation 

Wilson lines and their renormalisation properties are crucial ingredients in a variety of 
formal and phenomenological applications of quantum field theory. It is well known that 
the vacuum expectation value of a product of Wilson line operators, which we may denote 
generically by S, renormalises multiplicatively to all orders in perturbation theory [1-4]. 
As a consequence, such quantities S obey linear evolution equations of the form 



where T$ is the corresponding anomalous dimension. The solution of this equation has the 
form 



so correlators of Wilson lines generically have an exponential form. 

In an Abelian gauge theory, the anomalous dimension is simply a scalar function of 
the renormalized coupling and the kinematics and charges of the Wilson lines. In a non- 
Abelian gauge theory, things are more complicated due to the non-trivial colour structure. 
Specifically, for a correlator involving four or more eikonal lines joining at a hard interaction 
vertex, or equivalently a self-intersecting Wilson loop, there are several different possible 
contractions of the colour indices, corresponding to different colour flows at the hard vertex, 
as illustrated in a simple case in figure 1. Because gluons carry colour charge, these colour 
flows mix under evolution, implying that (1.1) is a matrix equation. The exponential 
solution for S in (1.2) is then defined through its Taylor expansion and the matrices in each 
term are ordered in correspondence with the scale A (this is indicated by the symbol V). 



Figure 1: The process qq — ¥ qq as a simple example of multiple colour flows: the two diagrams 
on the r.h.s. represent the two possible contractions of colour indices in the hard interaction, 
namely, c^ 7 = SijSki and c% = Si^Sji. Note that gluon exchange between the quark lines mixes 
between the two hard-scattering colour-flow components, leading to a matrix evolution equation in 
the two-dimensional colour-flow space defined by the two tensors C\ and ci. 

A complementary picture by which exponentiation can be understood is that of 'webs'. 
Considering the configuration of two eikonal lines meeting at a hard vertex (as exemplified 
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by the diagrams of figure 2), webs in an abelian theory are simply the set of connected 1 
diagrams as depicted in the figure. Multiple photon exchange diagrams between the Wilson 
lines, such as the ladder or crossed diagrams of figure 3, are generated by exponentiation [5]: 
they are fully reproduced by higher-order terms in the expansion of the exponential of a 
single exchange. Therefore, the exponent itself does not include any contributions from 
multiple-exchange diagrams, and its structure is very simple. 




Figure 2: Examples of connected diagrams in an abelian gauge theory; these diagrams exponenti- 
ate. The straight lines represent semi-infinite Wilson lines meeting with a cusp at the origin, where 
a charge-conserving hard interaction takes place. 




Figure 3: Examples of non-connected diagrams in an abelian gauge theory: these diagrams do not 
contribute to the exponent. The straight lines represent semi-infinite Wilson lines as in figure 2. 

In a non-abelian theory the situation is complicated: multiple interactions with a given 
Wilson line involve colour generators, and therefore do not commute. Nevertheless, it is 
known for the specific case of two eikonal lines (with a colour-singlet hard interaction) that 
exponentiation in terms of webs can indeed be generalised to the non- Abelian case [6-8]. 
That is, the relevant correlator of Wilson lines has the form 



wher J-(D) and C(D) denote, respectively, the kinematic dependence and the "Exponen- 
tiated Colour Factor" (ECF) of a given Feynman diagram D. This expression has the 

x Note that the eikonal lines themselves are excluded: subdiagrams with gluons that are attached to the 
same eikonal line(s) are not considered connected. 




D 



(1.3) 
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same schematic form as the abelian case of figure 2, but the contributing class of diagrams 
is larger, and has a more involved structure: this class comprises all "colour-connected" 
diagrams, namely the ones that cannot be partitioned by cutting only the eikonal lines. 
Two-loop examples for connected, colour-connected and non-colour-connected diagrams in 
a non-abelian theory are shown in figure 4. 

While the kinematic dependence of a given diagram, J-(D) in (1.3), is computed as 
usual using the eikonal Feynman rules, the corresponding colour factor C(D) differs from 
the original colour factor associated with that diagram, C(D). Gatheral [7] derived an 
iterative formula by which these ECFs C{D) can be computed, and characterized them 
as "maximally non-abelian". Soon thereafter Frenkel and Taylor [8] clarified the precise 
meaning of this term, defining a "web" as a "colour-connected" diagram. 




Figure 4: Examples of two-loop diagrams in a non-abelian gauge theory: (a) connected; (b) colour 
connected; and (c) non colour connected. The straight lines represent semi-infinite Wilson lines (here 
taken in the fundamental representation of the gauge group) meeting with a cusp at the origin, where 
a colour-conserving hard interaction takes place. The ladder diagram (c) does not contribute to the 
exponent, while diagrams (a) and (b) do. However, diagram (b) contributes with an exponentiated 
colour factor C(b) — —CaCf/2, instead of its original colour factor C(b) = Cf{Cf — Ca/2), whereas 
for diagram (a) the exponentiated colour factor equals the original one: C{a) — C{a). 

The conclusion of Ref. [6-8] is then that in the two-eikonal line case, the diagrams D 
that contribute to the exponent, that is the "webs", are the two-eikonal-line irreducible 
diagrams. We shall see that the definition of webs will have to be revised in the multi-leg 
case. 

Despite the fact that webs provide a direct graphical description of the exponent, they 
have so far been defined only in the particular case of a colour-singlet hard interaction 
with two eikonal lines 2 . The definitions, and indeed the formalism for computing the ECF 
of Refs. [7, 8] apply to a Wilson loop with any number of cusps or colour-singlet hard 
interactions, but they do not apply to self-intersecting Wilson loops which give rise to 
multiple colour flows. 

The purpose of the present study is to extend the concept of webs to multi-parton 
eikonal interactions, or equivalently to self-intersecting Wilson loops with multiple colour 

2 The three eikonal line case was analysed in Sec. 3.3 of Ref. [9]. Owing to the fact that for three lines 
(in the fundamental or adjoint representations) there is a single colour flow, this case ends up being a 
straightforward generalization of the two-eikonal line case, where only irreducible diagrams appear in the 
exponent. 



- 4 - 



flows. We will consider correlators of an arbitrary number of semi-infinite Wilson lines 
joining up at a single vertex (corresponding to a single hard interaction with an arbitrary 
colour exchange), and show that webs can be systematically defined, leading to exponen- 
tiation in the form of (1.3). Furthermore, we will present a general algorithm by which 
the ECF of any multi-loop (and multi-leg) diagram can be directly computed as a linear 
combination of ordinary colour factors of a definite set of diagrams, and present explicit 
results in several non-trivial cases. 

1.2 Motivation 

While our method is completely general, our motivation stems from concrete applications 
to scattering amplitudes. A major long-term goal is to have a complete understanding of 
long-distance singularities in multi-leg gauge-theory scattering amplitudes with arbitrary 
kinematics and colour representations. It is known that long-distance singularities expo- 
nentiate, where the exponent admits an all-order structure which is far simpler than the 
amplitude as a whole. The singularities therefore provide a window by which the all-order 
structure of perturbation theory can be explored. The general singularity structure is also 
universal amongst different gauge theories. In particular, a large class of singularities in 
any amplitude are controlled by the cusp anomalous dimension 7_ft-(a s ) [1, 10-14], asso- 
ciated with the renormalization of Wilson lines with a cusp. This object is by now well 
under control in QCD [15] , and even more so in the case of N = 4 super- Yang-Mills, where 
owing to integrability, a connection was established between the weak and strong coupling 
expansions [16-19]. 

Beyond the obvious field-theoretic interest in long-distance singularities, determining 
their structure is essential for precision collider physics. Indeed, the most difficult aspect 
of implementing cross-section computations with general kinematics beyond the tree level 
in QCD is the fact that infrared singularities only cancel after summing up real and virtual 
diagrams which are integrated over phase space with different numbers of partons. Know- 
ing the singularity structure is therefore a prior condition for any loop-level computation. 
Moreover, the singularity structure of gauge-theory amplitudes provides the key to resum- 
mation of large logarithms [20-57], which in many cases is essential for providing precise 
predictions for cross sections. So far resummation has mainly been applied to completely 
inclusive cross sections, where typically only two coloured partons participate in the hard 
interaction. LHC physics requires resummed calculations for less inclusive observables and 
more complex processes, where more partons participate in the hard interaction. Better 
understanding of long-distance singularities, and in particular of the correlation between 
momentum and colour flow, is essential for progress on this front. 

Over the past couple of years there has been much progress in understanding long- 
distance singularities of gauge-theory scattering amplitudes for both the massless [58-67] 
and massive [33, 43, 68-75] cases, by considering the factorization of fixed-angle ampli- 
tudes into hard, jet and soft components. The soft component, which captures all non- 
collinear long-distance singularities, is defined as a correlator of semi-infinite eikonal lines 
joined together at a single interaction vertex, as shown schematically in figure 5. The 
soft function S is, in general, a matrix in colour space [30,76-80], and it obeys an evo- 
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Figure 5: The soft function in the scattering amplitude involving five hard coloured partons. The 
hard interaction structure cannot be probed by soft (long wavelength) gluons, and therefore it 
reduces to a local vertex. Similarly, the partonic jets are not resolved and couple to the soft fields 
through eikonal vertices. 



lution equation of the form of eq. (1.1), leading to exponentiation in terms of the "soft 
anomalous dimension" matrix T$ (this object is also sometimes referred to as the "cross 
anomalous dimension" [4,81]). This anomalous dimension encodes all correlations between 
momentum and colour flow, which are only visible considering the full non-abelian the- 
ory at finite iV c . Concrete motivation to study this function beyond the two- loop level is 
provided by the fact that in the massless case its structure is strongly constrained [60,63] 
by both factorization, and a rescaling symmetry with respect to the Wilson-line velocities. 
These constraints point to a remarkable possibility, namely that all soft singularities in 
any multi-leg amplitude take the form of a sum over colour dipoles formed by any pair of 
hard coloured partons [60,62,63]. This very minimal structure, not involving any multi- 
parton correlations, is consistent with all explicit amplitude calculations done so far. Only 
very particular types of corrections, depending on conformally-invariant cross ratios, may 
occur at the three- loop order [60,63,64]. An additional type of multi-parton correlation 
may potentially show up starting from four loops, through higher Casimir terms. Further 
progress in studying the soft anomalous dimension for multi-leg scattering at three-loops 
and beyond requires direct understanding of non-abelian exponentiation, for which webs 
are the basic ingredient. 

A further motivation is provided by the prospect of resumming logarithms near kine- 
matic thresholds beyond the leading power in the momentum of emitted soft gluons, or 
equivalently, beyond the eikonal approximation [34,35,82-89]. A systematic classification 
of the structure of next-to-eikonal corrections has been given in [89] for the case of two 
hard partons coupled by a colour singlet interaction (the setup originally considered by 
Gatheral). It was shown that a subclass of next-to-eikonal terms formally exponentiates, 
and the relevant diagrams constitute a well-defined generalisation of eikonal webs. These 
results suggest that it is also possible to classify those next-to-eikonal contributions which 
exponentiate in the multiparton case, provided one can also classify eikonal webs. We will 
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see that this is indeed the case. Furthermore, webs are beneficial even in eikonal resumma- 
tion applications, where an anomalous dimension description already exists. Essentially, 
webs contain more information than anomalous dimensions (see [56]), namely finite parts 
as well as singular terms. Classification of these terms is crucial to their resummation, 
which can be phenomenologically important. 

Finally, another theoretical motivation to have a direct handle on exponentiation via 
webs is provided by the observation that exponentiation of gauge-theory amplitudes, at 
least in some cases, appears to extend beyond the singular terms, and apply to the am- 
plitude as a whole. An interesting case where complete exponentiation was conjectured to 
hold [90] is the class of maximally helicity violating (MHV) amplitudes in N = 4 super- 
Yang-Mills (SYM). While the very simple exponential form has later been shown to fail 
starting with the two-loop six-point function [91], N = 4 SYM perturbative amplitudes 
certainly admit an iterative structure also in their finite parts. A further step in reveal- 
ing the structure of these amplitudes was the discovery of a surprising duality between 
scattering amplitudes in momentum space and expectation values of Wilson loops taken 
in an auxiliary coordinate space [92-97]. Further, much progress was made on the same 
amplitude - Wilson- loop relation at strong coupling, see in particular Refs. [16,98,99]. 
Importantly, all of this profound theoretical progress has been limited to the planar limit 
(sec, however, [100]), and it would be very interesting to extend it beyond this limit, where 
the non-trivial aspects of colour flow reveal themselves. 

Our goal in this paper, as stated above, is to extend the notion of web exponentiation 
to the case of multiparton scattering, or equivalently to self-intersecting Wilson lines. To 
this end, we adopt the approach of [89], which uses path integral methods to relate the 
amplitude for soft gluon emission in a given process to a field theory for the soft gauge field. 
In this way soft gluon exponentiation can be easily formulated to all orders. For example, 
in the abelian case it is equivalent to the known exponentiation of connected diagrams 
in quantum field theory. Furthermore, the path integral method allows a straightforward 
generalisation of exponentiation properties to beyond the eikonal approximation. In the 
non- Abelian case, the field theory obtained in [89] for the soft gauge field was complicated 
by the non-commuting nature of the source vertices for soft gauge boson emission. It was 
still possible to classify which diagrams contribute in the exponent of the scattering ampli- 
tude, however, by using the replica trick, a combinatoric method often used in statistical 
physics applications (see e.g. [101]). The diagrams thus obtained were shown to be entirely 
equivalent to the webs of [7, 8] . These results suggest that a similar combination of methods 
will prove fruitful in analysing the multiparton case. We will see that this is indeed true, 
once one generalises the replica trick argument from the two parton case considered in [89] . 

1.3 Results and outline of the paper 

We may summarise our main results as follows. Firstly, we provide an algorithmic method 
for determining which graphs enter the exponent of eikonal multiparton scattering ampli- 
tudes, and for computing their exponentiated colour factors (ECFs), to arbitrary order 
in perturbation theory. Secondly, we derive a combinatoric formula for the ECF of any 
graph. Thirdly, we show that the general structure of the exponent of the eikonal amplitude 
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consists of sums of terms of the form 

W {ni ,n 2 ,...,n L ) = d W = E R DD' C(D') , (1.4) 

D D,D' 

where the indices D and D' label Feynman diagrams, C(D) and C(-D) denote, respectively, 
the conventional and the exponentiated colour factors of diagram D, and F(D) the kine- 
matic part. As is perhaps already clear from the structure of eq. (1.4), one finds closed sets 
of diagrams which mix with each other in the exponent. Each such set of diagrams can 
be labelled by the number of gluon attachments rij to each of the eikonal lines i = 1...L. 
Distinct diagrams D in the set differ only by the order of attachments of the gluons to 
each eikonal line. Associated with each such set (ni,n2, . . . , til) is a mixing matrix Rdd' 
describing how the exponentiated colour factors differ from the conventional ones. Com- 
binations of terms such as that shown in eq. (1.4) are the generalisation of webs from the 
two-eikonal lines to the multi-line case. The study of web structure is then equivalent to 
the study of the mixing matrices Rdd>- in particular, we observe that these matrices have 
the property 

J2 R DD'=0. (1.5) 
D' 

As discussed in section 5.4, this corresponds to the fact that the symmetric colour part 
of any diagram (i.e. that part of the colour factor that does not depend on the ordering 
of the soft gluon attachments to a given eikonal line) does not contribute to the exponent 
of the eikonal amplitude. This is the appropriate generalisation of the notion of antisym- 
metrisation of gluons previously noted at two loops. 

We also relate further remarkable properties of the matrices R to the cancellation of 
subdivergences in the exponent. We observe that R is an idempotent matrix, R? = R, and 
indeed it acts as a projection operator. Its eigenvalues are consisting only of and 1. These 
pick out linear combinations of kinematic factors, where the zero eigenvalue combinations 
correspond to those which are generated by the exponentiation of lower order webs, and 
thus do not contribute to the exponent. The eigenvalue one combinations build up in 
the exponent. We find that the corresponding linear combination of kinematic factors are 
special in that they encode intricate cancellation of subdivergences, allowing the webs to 
conform with the renormalization of the multi-eikonal vertex [102]. 

Finally, we argue that a subclass of soft gluon corrections exponentiates at next-to- 
eikonal order in multiparton scattering amplitudes. This generalises the previous results 
obtained for two parton scattering [89], and combines the path integral approach of that 
paper with the generalised replica trick argument developed in the present paper. 

The structure of the paper is as follows. In section 2, we summarise in more detail pre- 
viously known results on non-Abelian exponentiation, in particular the general singularity 
structure and the relationship between the soft anomalous dimension and web frameworks. 
In section 3 we set up the replica trick for multiparton scattering and establish a gen- 
eral algorithm for computing the ECF of a given diagram. We then use this in section 4 
to provide explicit combinatoric results for ECFs in terms of conventional colour factors, 
as well as an inverse relation. In section 5 we provide several demonstrations of how to 
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apply the new formalism; after showing how known results are reproduced, we analyse 
higher-order diagrams in perturbation theory and investigate the properties of the ECFs. 
Specifically, in section 5.3 and in Appendix C we compute all the non-trivial ECFs at three 
loops, and in section 5.4 we analyse an interesting four-loop class of diagrams and then 
discuss some general features of the mixing matrices Rdd 1 - In section 6 we examine the 
cancellation of sub divergences in the exponent. In section 7 we discuss the exponentia- 
tion of next-to-eikonal corrections in multiparton scattering. In section 8 we present our 
conclusions. 

2. Non-abelian exponentiation via evolution equations 

Exponentiation is a fundamental property of correlators of Wilson lines, not of scattering 
amplitudes. Nevertheless, it becomes relevant to the structure of amplitudes owing to the 
fact that amplitudes factorize, in such a way that their infrared singularities are captured 
by ultraviolet singularities of corresponding correlators of eikonal lines [10-14], a concept 
which also underlies the formulation of soft-collinear effective theory [103-105]. In this 
section we review the factorization of fixed-angle multiparton scattering amplitudes and 
then define and discuss the evolution of the corresponding "soft function" (or "eikonal 
amplitude") which will be the central object we analyse in the remainder of this paper. 

The qualitative difference between the general L-parton scattering case we consider 
here and the two-parton case where the singularity structure is fully understood, is the 
presence of multiple colour flows, as illustrated in figure 1 above. Let us therefore begin by 
introducing the notation for colour flow in the general case. The amplitude M describes 
the scattering of L hard (massless or massive) gauge particles with momenta pi, plus 
any number of colour-singlet particles, so it is characterized by L colour indices {a^}, 
i = 1,...L, belonging to arbitrary representations of the gauge group; we denote the 
corresponding generators by Tj. M can be decomposed into components by choosing a 
basis of independent colour tensors with the same index structure. We denote these tensors 
by (cx){ ai }, where K runs over all irreducible representations of the gauge group that can 
be constructed with the given indices {a«}. The decomposed amplitude in this colour- flow 
basis is 



where a s = a s {^) is the renormalized coupling in d = 4 — 2e dimensions. A simple example 
for such a decomposition is given in figure 1. 

From the soft physics perspective, hard partons are seen merely as classical sources of 
radiation, eikonal lines whose kinematics is fixed: they do not recoil upon emitting soft 
gluons. This is the essence of factorization. 

Having factorized the amplitude, the soft physics is described by the ultraviolet sin- 
gularities of operators composed of products of eikonal lines. These operators renormalize 
multiplicatively, involving the so-called "soft anomalous dimension" matrix (1.1). Upon 
solving the corresponding evolution equations one recovers the infrared singularities of 




(2.1) 



K 
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the amplitude. Much of what we know today about the all-order structure of infrared 
singularities is based on this formalism. 

In the following two subsections we briefly summarise the picture of factorization of 
scattering amplitudes in the massive and massless cases, and explain how it constrains the 
all-order structure of long-distance singularities, which are encoded in the soft anomalous 
dimension. In this discussion we use dimensional regularization throughout. Subsequently 
we shall address the relation with the complementary picture of exponentiation through 
webs, which we develop in the rest of this paper. 

2.1 The singularity structure in the massive case 

Next we turn to discuss the singularity structure of Mk using the basic tool of factorization 
into hard, jet and soft components [59-63,76,78,106-108,108-110]. For simplicity consider 
first the case where all the hard partons are massive. In this case long-distance singularities 
only arise from the exchange of soft gluons (gluons with vanishing energy), and power 
counting shows that there is at most one such singularity per loop, i.e. at 0{a™) the 
renormalized amplitude has poles C(l/e fc ), with k = l...n. At each order, only the 
0(\/e l ) pole is "new", namely computing it requires an C(a") calculation, whereas all the 
higher poles can be determined from lower orders using evolution equations. 



Figure 6: Factorization into soft and hard functions in the example of four parton scattering with 
momenta The four lines represent Wilson lines along the directions pi in the colour representation 
corresponding to the respective hard parton. The labels bi, for i = 1 ... 4, on the external lines are 
the colour indices of the cikonal amplitude, while the corresponding labels a,; arc the colour indices 
of the hard subprocess H at the origin. 

If we assume that all kinematic invariants pi ■ pj , including the masses pf = mf , are 
taken simultaneously large compared to Aq CD , soft gluons completely decouple from the 
hard interaction: soft singularities are entirely independent of the overall hard scale. We 
may then factorise the amplitudes into a soft function 5, a matrix in the colour-flow space 
defined above, acting on the hard subprocess H, which similarly to M is a vector in this 
space. The essence of this factorization is the fact that all the singularities of the amplitude 
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are captured by S, while H is finite. This separation is useful because S has a much simpler 
structure as compared to the full amplitude; for example, it does not depend on the spin 
of the scattered particles. 

As illustrated in figure 6, the factorized amplitude takes the form: 



Mj (pi,a s ,e) = 



K 




1717 



2 > a sj e 



(2.2) 



where the soft function is defined as the vacuum expectation value of the following product 
of Wilson lines, 



(cj) {ak} S JK 



,a s (/x 2 




1^ 



) (c K ) {bk} . (2.3) 



The colour index structure of the Wilson lines is illustrated in figure 6. In (2.3) we projected 
the Wilson line correlator onto the colour- flow basis defined in terms of the tensors {ck}- 
Each Wilson line is defined as usual, as a path-ordered exponential: 



^Pexp 


ig s / dzfAp 


J = (-Pexp 




L Jci 





/ 

J 



\g s / dt^ ■ A(tft 



(2.4) 



going along the classical trajectory Ci of the hard parton I, that is starting at the origin 
and going in a straight line, with 4- velocity ft proportional to its momentum p\. The first 
argument of S indicates the dependence of this function on the kinematics: given L legs, 
it depends on L(L — l)/2 Lorentz invariants. The ratio (ft • Pj) 2 / ((3f f3j) is the hyperbolic 
cosine of the cusp angle formed by the two Wilson lines ft and ft in Minkowski space- 
time. The dependence on this ratio - rather than on ft • ft directly — is dictated by the 
invariance of the eikonal Feynman rules with respect to rescaling the velocities. 

The soft function S is defined as an operator composed of Wilson lines, and therefore 
it is multiplicatively renormalizable [3,4]. However, being a matrix in colour space its 
renormalization involves mixing [4], so the evolution equation takes a matrix form: 




,a*(M 2 )> e 



Sji 



(ft -ft) 2 
Pi 



,a s (^ 2 ),e T 



IK 



(ft -ft) 2 



,«s(^ 2 ) 



(2.5) 

Here Tf K is the "soft anomalous dimension" (it is called the "cross anomalous dimension" 
in Ref. [4,81]). It is a finite function of the kinematics and the d = 4 — 2e dimensional 
coupling constant. As shown in [111], evolution equations of this type can be solved directly 
in dimensional regularization: the singularities in S are generated upon integrating the 
evolution equation (2.5) from /i 2 = 0. Thus, as anticipated, at each order there is one new 
coefficient to determine, the C(a") term in V s , which would generate the 0(l/e v ) pole in 
S at this order in perturbation theory. 

Very recently the soft anomalous dimension in the massive case, T s , has been explicitly 
computed to two-loop order [43,68-72,74,75], where a direct correlation between the colour 
and kinematic degrees of freedom for three partons first appeared. 
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2.2 The singularity structure in the massless case 

Let us turn now to discuss the singularity structure in the massless case. Here an ad- 
ditional source of long-distance singularities is present, namely gluons that are collinear 
with one of the hard massless partons. As a consequence, there are two extra singularities 
per loop: that is, at a given order C(a"), the amplitude contains 0(l/e k ) singularities, 
where k = l...(2n). Here the real power of exponentiation becomes apparent: in the 
exponent for each extra power of a s there is just one extra singularity. At one loop the 
exponent has an 0{l/e 2 ) singularity from overlapping soft and collinear divergences, as 
well as subleading 0{l/e l ) singularities from either soft or collinear divergences; at order 
it contains 0(l/e k ) where k = 1 . . . (n + 1). The absence of 0(l/e k ) singularities for 
n + 1 < k < 2n in the exponent corresponds to the fact that multiple exchange diagrams 
(e.g. ones where each gluon produces a soft as well as a collinear divergence) are not part 
of the exponent, but are instead reproduced by expanding it. Similarly to the massive 
case, the only "new" singularity at any given order, which cannot be determined using the 
previous orders (knowing the cusp anomalous dimension and the /3 function coefficients) is 
the 0{l/e l ) pole. As in the massive case, this structure is dictated by evolution equations. 

To present the factorized amplitude in the massless case, let us assume for simplicity 
that all the hard partons involved are massless, p 2 = for any leg I = 1 . . . L, and take 
the remaining kinematic invariants pi • pj to be simultaneously large compared to Aq CD 
(the fixed-angle scattering limit). Long-distance singularities can then be shown to emerge 
from separate collinear and soft regions; these decouple from the hard interaction as well 
as from each other (the latter step is rather subtle, see for example [112]), leading to the 
following factorised structure [59,60,64-66,78]: 

/ 2p^ ■ p ■ (2pi • rij) 2 \ 
Mj (pi,a s (p 2 ),e) = 22 SjK ' ft' a s(/i 2 ),f) H K y ^ 2 3 , -^-j ,a s (p 2 ),ej 

jj&^.a.^.e) (2-6) 



n 



n 2 p 2 



i=1 ^i( m J i)2 ^ 2 ),e 

where the soft function is defined as in the massive case, (2.3), replacing each hard parton 
by a corresponding Wilson line: 



(cj) {ak} S JK (ft ■ (3 V a s (p 2 ), e)^J2 ( 

{b k } \ 

The evolution equation corresponding to S is: 



n [$A(°o,o) afcfefc 
i=i 



(c K ) {bk} . (2.7) 



^Sjk (ft • Pj, a s (p 2 ), e) = - Sjj (ft • ft, a s (p 2 ), e) Fj K (ft • ft, a s (p 2 ), e) . (2.8) 

There are two related aspects by which this equation differs from the one corresponding to 
the massive case (2.5): the first is that the anomalous dimension r 5 itself is now singular 
for e — > owing to the collinear singularities, and the second is that this function, and thus 
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S itself, depends on ■ (3j — not on the scale invariant ratio we encountered in the massive 
case — which implies violation of the rescaling symmetry with respect to the velocities 
The relation between these two aspects was explained in [60] : the singularities appearing in 
T s and the violation of the rescaling symmetry are both consequences of the cusp anomaly 
for massless Wilson lines, and they are both controlled by the cusp anomalous dimension 
1k(&s) to all orders. 

In order to define the partonic jet function Jj in (2.6) we introduce an auxiliary vector m 
(such that nf / 0) [113] and form the following gauge-invariant transition amplitude [59,60]: 

u(p) J ( {2P 2 n f ,a s (fi 2 ),e) = (p | ^(0) $ n (0, -oo) |0> . (2.9) 

The overlapping soft and collinear divergences to jet i are contained in both S and Jj. 
Therefore, to avoid the double counting of this overlap region we have divided each jet J, L 
in (2.6) by its eikonal part, Ji, which is already contained in S. The eikonal jet is defined 
by: 

J ( 2(/3 n ' 2 n) ^s{n 2 ), e) = <0| <Moo, 0) $ n (0, -oo) |0) . (2.10) 
It follows then that the "reduced soft function" [59, 60] 

c ( / 2\ \ Sjk (Pi ■ Pj, at s (fJ> 2 ), e) 

Sjk (pij,a s {fi ),e) = — '- (2.11) 

^ / o ( a . ^ . \ 2 



^ ( 2(p i .n l f 2 \ 



i=i v * 

is free of collinear poles, and contains only infrared singularities originating from soft gluon 
radiation at large angles with respect to all external legs. This function admits a similar 
evolution equation to (2.8), 

Vj^S-JK {Pij,oi s (n 2 ),e) = - Sji (pij,a s (fi 2 ),e) Tf K (p i: ,-, a s (fi 2 )) . (2.12) 

but with an important difference: the anomalous dimension matrix V s is now finite, cor- 
responding to the fact that S has no collinear singularities, just soft ones. 

In S the rescaling symmetry is restored, and thus it depends exclusively on the following 
kinematic variables 3 : 



(-fr-Pj) 2 IPi-Pjfe' 2 ^ 



" 2(p t -n t ) 2 2(P J -n J ) 2 2(ft • n-) 2 2{P j ■ n 3 ) 2 ' {2A ' i} 



rcf n 2 n 2 n 2 



which are inherently scale invariant. Ref. [60] constrained V s in (2.8) by enforcing the 
cancellation of rescaling violation between the numerator and denominator in (2.11). This 
led to a set of L all-order matrix equations [60]: 

E diip—) r ^ {Pij > as) = \ 7 * {as) SlK ' v * ' (2 - 14) 



3 The phases Xij are defined by f3t ■ /3j = — ■ /3j\e 17rAij , where Xij = 1 if i and j are both initial-state 
partons, or both final-state partons, and Xij — otherwise. 
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which constrain the kinematic and colour dependence of the soft anomalous dimension of 
any massless scattering amplitude. An equivalent set of constraints was derived in [63] 
using the soft-collinear effective theory formulation for the problem. 



As mentioned in the introduction, the simplest solution to eqs. (2.14) is the sum-over- 
dipoles formula [60,62,63]: 



where and 5$ are eikonal anomalous dimensions, depending exclusively on the d- 
dimensional coupling (see [60]). The first term in (2.15) correlates directly the colour 
and kinematic degrees of freedom of pairs of hard partons, while the second term is colour 
diagonal, and can be absorbed into the jets. An important aspect of this solution is that 
the entire matrix structure in colour space is governed by the cusp anomalous dimension 
to all orders in perturbation theory. 

Eq. (2.15) is indeed a solution to (2.14) if and only if the cusp anomalous dimension 
admits Casimir scaling, i.e. j^{a s ) = Tj • Tj (a s ), where (««) does not depend 
on the representation of the parton i. If the latter condition is violated at some order in 
perturbation theory - and, as discussed in [60,64], this may happen starting at four loops 4 
- there will be additional contributions to V s involving quartic and higher Casimirs. 

In addition, solutions of the homogeneous equations corresponding to (2.14) may con- 
tribute to V s starting at three-loops [60] . These are functions of conformally- invariant cross 
ratios of the form PijPkif(PikPjl), which simultaneously correlate the colour and kinematic 
variables of four or more partons. A dedicated analysis of such corrections was performed 
in [64] taking into account a range of additional constraints, including Bose symmetry, 
collinear limits [63] and considerations of the degree of transcendent ality. This analysis 
showed that although the class of possible functions is strongly constrained, such correc- 
tions cannot be altogether excluded even at the three loop order. Moreover, at three-loops 
there is a unique logarithmic function of conformally-invariant cross ratios that fulfils all 
constraints, possibly providing a first genuine multi-parton correlation in the massless case. 
An explicit 3-loop 4-leg calculation would be necessary to determine whether such a cor- 
rection to (2.15) is indeed present. Such a calculation would greatly benefit from methods 
for efficiently classifying which diagrams contribute to the exponent, which is the subject 
of this paper. 

2.3 Anomalous dimensions and webs 

In the previous subsections we have briefly summarised the conventional approach to non- 
abelian exponentiation in terms of evolution equations. The first conclusion from this 
discussion is that in both the massive and massless cases, the physics of non-abelian expo- 
nentiation is encapsulated by the correlator of semi-infinite Wilson lines, defining the soft 
functions in (2.3) and (2.7) respectively. It is this correlator that will be the subject of our 
diagrammatic analysis in what follows. 

4 See however [63] where evidence is provided against such four-loop contributions. 



r 5 {Pij,a 9 ) 



1 



^ (<*«) E E Mfti) Ti ■ Tj + X - Sg(a a ) ^ T, • T, , (2.15) 



dip. 



8 
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Studying soft anomalous dimensions we saw that the kinematic dependence of such 
Wilson line correlators is tightly linked to their structure in colour space, and that this 
interesting physics is only visible in the non-abelian theory at finite N c . 

Finally, we learnt that non-collinear soft singularities are always described by finite 
anomalous dimension matrices, through eqs. (2.5) and (2.12) in the massive and massless 
cases, respectively. This corresponds to the simple but crucial fact mentioned above, that at 
each order in perturbation theory only the 0(l/e l ) pole in the exponent encodes genuinely 
new information, which cannot be deduced from lower orders. In terms of Feynman dia- 
grams this observation is highly non-trivial: after summing up the diagrams corresponding 
to the exponent of (2.3) or (2.7), and accounting for the renormalization of the coupling 
(and in the massless case, for the eikonal jets) all remining singularities are associated 
exclusively with the renormalization of the multi-eikonal vertex. The latter comprises at 
each order, of just one new counter term removing a single ultraviolet divergence. Subdi- 
vergences must therefore conspire to cancel. 

In the two-parton case (colour singlet hard interaction) of [6-8], the situation is simple: 
diagrams that contribute to the exponent just do not have any subdivergence. So upon 
considering the massive case, where no collinear singularities are present, and accounting 
for the renormalization of the coupling, the diagrams that appear in the exponent indi- 
vidually have a single pole each, while others, such as the ladder diagrams, which have 
subdivergences, do not contribute to the exponent. For example, in the massive case, at 
two loops, the diagrams in figure 4 (a) and (b) have no subdivergences, as can be seen 
from the fact that shrinking one of the gluon loops towards the cusp drags in the other 
attachments, leading to a single overall ultraviolet divergence, while in (c) there is a sub- 
divergence (the internal gluon can shrink to the cusp without affecting the external one) 
and indeed this diagram does not appear in the exponent. 

In the multi-parton scattering case subdivergences are obviously present in individual 
diagrams that contribute to the exponent. This implies that intricate cancellation must 
take place between diagrams. This, as we shall see, goes hand in hand with the need to 
generalise the concept of webs in the multi-leg case from individual diagrams to particular 
combinations of diagrams that mix with one another as described by (1.4). We shall 
return to this fundamental issue in section 6, after learning the properties of multi-leg 
webs. As discussed in the introduction, our analysis will make use of the replica trick, and 
we introduce this idea in the next section. 

3. Non-abelian exponentiation via webs using the replica trick 

In this section, we introduce the replica trick for multiparton scattering, which we will 
later use to study the structure of the exponent at higher orders. Our method is inspired 
by similar methods in statistical physics problems (see e.g. [101]), and was already used 
in [89] to derive the structure of eikonal and next-to-eikonal webs in two parton scattering. 
The argument presented in that paper is not immediately generalisable to the L parton 
case. Thus, here we reformulate the replica trick argument directly for the general L-leg 
case, such that two (and three) parton scattering emerge as special cases. For the moment 
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we focus on the eikonal approximation, and postpone the discussion of next-to-eikonal 
corrections to section 7. 

This section has two parts. In the first we present the formalism, based on the replica 
trick, for solving the combinatorial problem of non-abelian exponentiation in the multi-leg 
case. In the second we present an algorithm for computing ECFs based on this formalism. 

3.1 The replica trick formalism for non-abelian exponentiation 

Our starting point is to consider a hard scattering process H(x±, . . . ,xl) which produces L 
hard outgoing particles with initial positions {xi} and final momenta {pi}. In the eikonal 
approximation, one may approximate the spacetime path Zk(t) of the k th hard external 
particle by its classical trajectory 

z k (t) = x k + /3 k t, (3.1) 

where (3 k is the velocity factor already introduced in the previous section, and the initial 
positions all coincide at the hard interaction vertex, so that x k = 0. Then the scattering 
amplitude for this process, dressed by any number of virtual soft gluons can be written as 
follows: 

M bl ... bL (pi, ...,Pl)= [ [T>A*\ H ai ... aL (0, . . . ,0)e iS ^ ]J (Vexp \ig s f dtp k ■ A S V\ . 

J k V L J J / a k b k 

(3.2) 

This expression contains a path integral over the soft gauge field As (matrix-valued in 
colour space), whose action is denoted by S(Ag ). Each external line contributes a path 
ordered exponential, a Wilson line operator in the direction f3 k , where a k and b k represent 
the initial and final colour indices on line k, as illustrated in figure 6. 

Equation (3.2) generalises the expression given in [89] for the two-parton amplitude 
to any number of legs. It is a reformulation of a well-known result: outgoing particles 
emitting soft radiation in physical scattering processes can be treated as Wilson lines (see 
e.g. [23]). Using the same path-integral method Ref. [89] also considered what happens 
beyond the eikonal approximation, which we return to in section 7. 

One may take the hard interaction outside the path integral and rewrite eq. (3.2) as 



Mb 1 ...b L (Pl, ■ ■ - ,Pl) = Ha x ...a L Za x ...a L ,b 1 ...b L , (3.3) 



where 



Z ai ...a L , bl :.b L = j [^]e i5 ^n(Vexp ig s j dt (3 k • A 



(3-4) 

k \ l J J / a k b k 

may be recognised as the generating functional of a quantum field theory for the soft gauge 
field. The Wilson line factors act as source terms for the soft gauge field, and generate 
eikonal Feynman rules for gluon emission from the emitting particle legs. Carrying out the 
path integral thus generates all possible soft gluon subdiagrams which span the external 
lines. Using the Wilson line notation of eq. (2.4) for each line I = 1,2, ... L one may write 
eq. (3.4) more compactly as: 



s ^ [$« ® $( 2 ) ® . . . ® , (3.5) 
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where the tensor product symbol reminds us that the Wilson lines operate on different 
parton legs, involving the respective colour indices in the corresponding representations 5 . 

To address soft-gluon exponentiation, we essentially have to compute InZ. Here the 
replica trick comes in handy, which proceeds as follows. One considers N identical copies 6 
or replicas of the soft gauge boson field, which all share the same gauge group, say SU(iV c ). 
Importantly, the replicas do not interact with each other, so that the soft gauge field actions 
corresponding to different replicas <S[-AJJ combine additively, 



N 



(3.6) 



where A^ is the gauge field associated with replica number i (we have dropped the sub- 
script s denoting softness). Note that if there is matter in the vacuum, as in QCD, it is 
replicated as well. Thus, different replicas cannot communicate through matter (e.g. via 
gluons coupling to a fermion bubble), so they all have the same standard renormalization 
properties. 

By analogy with eq. (3.5), the generating functional for this theory is given by 



J [VAl]...[VA^] jXM] ($( 1 )$( 1 )...$«) (c|,( 2 )4 2 ) 



(2)> 
N i 



(3.7) 



where <J>^ is the Wilson line factor associated with parton I (where I = 1 . . . L) and replica 



number i (where i = 1. . .N). Note that only provides a source for the field A^ 
the particular replica i, however, owing to the non-abelian nature of the theory it does not 
commute with other $^ (where j ^ i). 

On the left-hand side of (3.7), we have recognised that the generating functional for 
the replicated theory constitutes that of eq. (3.5) raised to the N th power. Each external 
line now carries a product of Wilson lines, which has the form 



of 



aifei 



V exp 



J / aiC2 



x V exp 



1.9. 



Pexp 



19, 



C2C3 



(3.8) 



In order to derive the Feynman rule for emission of a soft gluon from parton I, the product 
of path-ordered exponentials must be combined into a single path-ordered exponential. 



5 Note that it is straightforward to project Z onto a specific colour-flow basis (cj) as done in (2.3) and 
(2.7), obtaining: 

Zu = f [VA*] e iS[A = ] ® $ (2) ® • • • ® $ (L) 



where we have associated the indices I and J with (ai . . . a^) and (61 . . . bi), respectively. 

6 The number of replicas N is a new integer parameter (it should not be confused with the number of 
colours N c ). This parameter is only introduced as a tool to handle the combinatorial problem of non-abelian 
exponentiation, and it need not be thought of as having a physical meaning. 
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This can be accomplished (in analogy with the two eikonal line case of [89] ) by introducing 
an operator 1Z that orders the colour generators associated with emissions along external 
line I according to their replica numbers. That is, one has 



= 1ZV exp 



i=l 



(3.9) 



a-lh 



where the sum in the exponent is over replica numbers 7 . Given that we have chosen all 
momenta to be outgoing, the operator 1Z reorders the colour matrices on a given line 
such that the replica number increases as one moves away from the hard interaction. 
Equation (3.7) thus becomes 



Z N = j [2MJ] . . . e^^yfcjpexp ig^Jdt^^ 



V exp 



N 

i=i j 



(3.10) 



where H reorders the colour generators along all the lines. 

Carrying out the path integrals over the soft gauge fields in (3.10) generates all soft 
gluon diagrams — involving all possible replica numbers — which connect the external 
lines. Since the actions -SLA*] are all the same as in the original (unreplicated) theory, the 
kinematic part of each diagram, J-~(D), is the same as one would obtain in the unreplicated 
theory for the same configuration of gluons (i.e. the operator 1Z acts only on the colour 
generators of gluon fields, and leaves their spacetime dependence unchanged). However, 
the colour factor of each generated diagram is not the same as would be obtained in the 
unreplicated theory (C(D)). Rather, the colour factor of a given graph in the replicated 
theory, which we may denote by Cn(D), is that which results after reordering the replicas 
according to the 1Z operator, and it explicitly depends on the number of replicas N. We 
will see explicit examples of this in section 5. 

Next, one may perform a Taylor expansion in N to obtain 



Z N = 1 + NlnZ + 0(N 2 



We now have 



]nZ = ^T(D) C(D), 

D 



(3.11) 



(3.12) 



where the sum on the right-hand side is over all diagrams D whose colour factors are 
0(N). Each diagram contributes the kinematic part T(D) as in the original theory, times 
a modified colour factor C(D) which is defined as the coefficient of N 1 in the Taylor 



7 Note that in the absence of 1Z on the r.h.s. the expansion of the path-ordered exponential would 
yield terms with all possible orderings of the colour matrices corresponding to the different replicas, not 
reproducing the l.h.s. where this order is uniquely fixed. Thus the 1Z operation is crucial for this equality 
to hold. 
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expansion of the colour factor in the replicated theory (note that C{D) itself carries no iV 
dependence). The sum in (3.12) is over all diagrams that have a component that is linear 
in the replica number N. Indeed, some diagrams do not have any 0(N 1 ) terms in their 
Taylor expansion, and therefore they do not contribute to the exponent. We shall see such 
examples in section 5. 

Given that eq. (3.12) is manifestly independent of N, the final result for the generating 
functional Z in the original (unreplicated) theory is therefore, 



which is manifestly an exponential form 8 . 

3.2 A replica-trick based algorithm for computing exponentiated colour factors 

The formal manipulations leading to (3.13) are useful because they uniquely specify the 
exponent: they directly translate into a general algorithm by which the ECF associated 
with any given diagram D can be computed. We may summarise this algorithm by the 
following sequence of steps: 

1. Identify distinct connected pieces in the diagram (see footnote 1 in section 1.1). We 
assume there are n c such pieces, fi, with i = 1 . . . n c . If n c = 1 the diagram as a 
whole is connected and C(D) = C(D). For n c > 2, assign each connected piece fi in 
the diagram a replica variable rj, which can take (integer) values from 1 to AT. 

2. Compute the colour factor Cn(D) in the replicated theory. To this end one needs to 
sum over all possible assignments of the variables rj in the range 1 to N. In practice, 
all that matters, based on the definition of the ordering operator 1Z, is the hierarchy 
of r.j amongst a set of gluon attachments to a given leg. Specifically, colour generators 
associated with two gluon attachments to a given leg will be ordered such that the 
generator with the higher replica number is further from the hard interaction. Given 
two attachments i and j to the same leg I, the 1Z operation 

• does nothing if rj = rj, in which case the colour factor will be determined by 
the ordering of the gluons in the original diagram D. 

• will order them as TjTj if rj < rj and as TjTj if rj < rj, independently of the 
order of the two attachments in the original diagram D. 

A convenient way to compute Cn(D) is therefore: 

8 A comment is in order regarding to the nature of the replica trick as used here. In statistical physics 
applications, the replica number N typically becomes a complex parameter. Recovering the unreplicated 
from the replicated theory often becomes ambiguous, due to subtleties involving analytic continuation. 
Here there are no such problems, the colour factor in the replicated theory is a polynomial in N, and it is 
straightforward to take its ©(JV 1 ) coefficient. 




(3.13) 
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I. Tabulate all possible hierarchies which may arise when each of the variables 
ri takes values in the range 1 to JV. A hierarchy determines whether r\ > rj, 
ri = rj or ri < rj for any two connected pieces i and j. For a given hierarchy h, 
we denote the number of different rj's (different active replicas) by n(h). This 
number can vary from 1 to the number of connected pieces n c . An example is 
provided in table 1 for the n c = 3 case (note that for simplicity the variables 
rj,7j and are denoted there by i, j and k, respectively). 

II. For each hierarchy in the table compute its multiplicity, namely a combinatorial 
factor counting how many times this hierarchy arises when all rj go over the 
range 1 . . . N. The multiplicity is equal to the number of ways of choosing n(h) 
different replica numbers from N (associating them to the n(h) different ?Vs), 
thus is given by the combination function: 

III. For each hierarchy h, determine the colour factor of the diagram after applying 
the ordering operation TZ to the generators on all legs (TZ is only relevant for 
legs to which there are two or more attachments). We denote this colour factor 
by TZ [C(D) \h], to be read as "the ^-ordered colour factor of diagram D, given 
the hierarchy h" . 

IV. Sum over all possible hierarchies h with their respective multiplicities M^(h), 
to obtain the colour factor of the given diagram D in the replicated theory as a 
function of N: 

C N {D) = Y, M N (h) TZ [C(D) \h] (3.15) 
h 

3. Expand Cjv(-D) in powers of N and extract the coefficient of iV 1 , which is the ECF 
C(D). 

One can go further and recast this algorithm as a formula for the ECF. The crucial 
observation here is that in our expression for Cn(D) only the multiplicity factor Mj^{h) 
depends on N, and this dependence is through the simple combination function in (3.14). 
We can therefore extract the iV 1 coefficient, obtaining: 

_^ (-l)n(h)-l 

C ^ = E i {h) *> [°( D ) N • ( 3 - 16 ) 

We shall continue to develop this in the next section, where we provide a more explicit 
formula for the ECF. 

Implementing the above algorithm requires a notation that identifies connected pieces 
and simultaneously specifies the order of attachments of gluons to each of the legs (eikonal 
lines) in an arbitrary diagram. Considering a set of diagrams 

{D} = (ni,n 2 ,n 3 , ■ ■ ■ ,n L ) 
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with L legs and ni gluon attachments to a given leg I, we use the following notation for a 
particular diagram within this set: 



D 



Each of the square brackets corresponds to a leg (there are L such terms in D). Each 
of the entries within a given bracket [sf \ ■ ■ ■ Sn}] corresponding to leg /, is associated 
with an individual gluon attachment, where the order of the list indicates the order of 
gluon attachments to this leg: the list is ordered from the outside inwards toward the hard 
vertex (cusp). The variables sf^ themselves are assigned values from 1 to n c (the number 
of connected pieces in D) associating each of the gluons with a particular connected piece 
fi in D. This facilitates the assignment of replica numbers. Note that this notation 
does not specify the way by which the gluons are connected internally (which may involve 
additional loops). These details are irrelevant to the combinatorial problem of non-abelian 
exponentiation. The notation (3.17) is used in section 5 where we consider explicit examples 
- see figures 14 through 19, and also Appendix C. 



4. General formula for Exponentiated Colour Factors 

In the previous section we have set up the replica trick formalism for multiparton eikonal 
scattering. This allowed us to solve the combinatorial problem of non-abelian exponentia- 
tion, constructing directly the exponent of such amplitudes in terms of individual diagrams. 
We have seen that the exponent contains those diagrams which, in the replicated theory, 
have a term that is linear in the number of replicas N. The contribution of each such 
diagram D to the exponent is determined by an ECF C(D). We have seen that the replica 
trick formalism translates into a general algorithm for computing these ECFs. We will now 
show that it also leads to an explicit formula for the ECF of a given graph D in terms of 
a particular linear combination of the products of ordinary colour factors of subgraphs in 
D. This is the multiparton generalisation of the similar result given in [89] for the case 
of two eikonal lines. Furthermore, in section 4.2 we present an inverted result, expressing 
the conventional colour factor C(D) of graph D in terms of ECFs. The latter result is of 
conceptual interest as the multiparton generalisation of Gatheral's formula for webs in two 
parton scattering [7] . 

4.1 Exponentiated Colour Factors in terms of conventional ones 

In order to proceed, let us first introduce the notion of decompositions 9 of a diagram D. 
A decomposition is a complete set of subgraphs of D - not necessarily connected - each 
containing a single replica number. A given decomposition P contains a number of replicas 
n(P), which in general is not sufficient to uniquely define the decomposition. A two-loop 
example is given in figure 7. Each square bracket represents one decomposition, and there 
are two in total with n{P) = 1 and n(P) = 2, respectively. In this simple case, there is only 
one decomposition for each choice of the number of replicas. Another example is shown in 
figure 8. Here there are 5 decompositions in total. 
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(a) (b) 

Figure 7: Decompositions of an example two loop diagram. The two decompositions have (a) 
n(P) = 1; (b) n{P) — 2. Different colours represent different replica numbers. 




Figure 8: Decompositions of an example three loop diagram, with (a) n(P) = 1; (b) n(P) = 2; 
(c) n(P) = 3. Different colours represent different replica numbers. 

Each decomposition contains a number of subgraphs g\ , 52 , ■ • ■ g n ( P) > each of them 
having a different replica number 10 . We now wish to express the ECF C(D) in terms 
of products of ordinary colour factors of these subgraphs C(gi). To this end we need to 
determine the colour factor in the replicated theory, Cn(D). We shall see that owing to the 
replica ordering operator TZ, the latter contains a sum of products of C(gi) in all possible 
orders. 

The simplest way to determine this sum is to use eq. (3.15) and split the sum over 
hierarchies h in this equation into a sum over decompositions P, where for each decomposi- 

9 Note that these were called partitions in [89]. 

10 The replica numbers of gi and gj (i 7^ j) in a given decomposition P cannot be equal. In this formalism 
the case where the two corresponding subgraphs are assigned the same replica number is taken into account 
in some other decomposition P' in which the two belong to one subgraph g^. 
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tion there is an internal sum over possible hierarchies. Crucially, this internal sum becomes 
trivial, since given the decomposition P all the possible hierarchies have n(h) = n(P), and 
therefore they have identical multiplicities Mn{K), which we denote Mn(P). Given P, for 
any particular h the 7£-ordered colour factor would be a particular product 

K [C(D) \h] = C(g mi ) . . . C(g n{hUp) ), (4.1) 

where ir(h) is the permutation of (1,2, . . . ,n(P)) corresponding to the hierarchy h. But 
since all M^{h) given P are the same, all the permutations contribute symmetrically to 
Cn(D). One therefore obtains: 

C n (D) = J2Mn(P)J2 C(g n )...C(g nn{p) ), (4.2) 

P 7T 

where the internal sum goes over all permutations of (1,2,..., n(P)) and Mjy(-P) is given 
by the combination function in (3.14) with n(h) = n{P). 

We can now proceed, as in eq. (3.16), to extract the 0(N 1 ) part of eq. (4.2). We 
obtain: 

(_-\\n(P)-l 

C(D) = £ 1 > £ C(g n ) . . . C(g nn(p) ), (4.3) 

P \ ' TT 

which thus represents the ECF of graph D, written explicitly in terms of ordinary colour 
factors of its subgraphs. A couple of examples for the application of eq. (4.3) are provided 
in appendix A; the same examples will be analysed in section 5 by applying the replica 
trick algorithm according to section 3.2, as well as by direct exponentiation. 

It is interesting to note a further simplification of this result when the colour factors 
C(X) of subgraphs commute with each other. When this is the case, one may write 

£ C(g n ) . . . C(g^ (p) ) = n(P)\ \{ C(g), (4.4) 

tt geP 



so that eq. (4.2) becomes 



C N (D) = ^2 N P n(P) l[C(g), (4.5) 

P geP 



where 

N\ 

NPn ^ = (N-n(P))\ (46) 
is the permutation function. Expanding in powers of N, 

N P n(P) = (-l)^)" 1 (n(P) - 1)! N + 0(N 2 ) (4.7) 

we extract the ECF for the case of subgraphs with commuting colour factors: 

C(D) = J^-l)^- 1 (n(P) - 1)! [J C(g) . (4.8) 

P g&P 

This is precisely the result that was given in [89] for the two eikonal line case, when the 
colour factors do indeed commute. 
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4.2 Conventional colour factors in terms of exponentiated ones 

The above formula expresses the exponentiated colour factors in terms of the conventional 
colour factors. It is also instructive to write an inverted formula which expresses things the 
other way around. Such a formula is the generalisation of Gatheral's original definition of 
exponentiated colour factors [7] to the multileg case. 

The argument here is similar to that presented in [89] for two eikonal lines, and runs 
as follows. One again uses the notion of decompositions of a diagram D into subdiagrams, 
as described above. A given decomposition is labelled by a set of numbers {mn}, where 
ran counts how many times a given subdiagram H appears in the decomposition. The 
index H runs over all possible subdiagrams. For example, the decomposition of figure 8(a) 
has nisa = 1 and all other m# = 0; figure 8(c) has m\ a = ran, = m\ c = 1 and all other 
mg = (n.b. we have used labels introduced in figures 13, 15, 16 and 17). One may write 
the exponentiated amplitude as 

exp ( £ F{H) C(H) \ = £ I ( £ F{H) C(H)) , (4.9) 

I H J n=0 V H / 

where as usual J~(H) represents the kinematic part of diagram H. The expansion of the 
exponential on the right-hand side contains a product of sums. One may rewrite this as 
a sum over all possible products of subdiagrams by using the notation for decompositions 
introduced above: 

(^T(H)C(H)) = £ ( J] HH)m " ) [c(ff 1 ) mi C(g 2 ) m2 . . . + perms] . (4.10) 

V H J { m „} V H mH ' / 

Some explanatory remarks are in order. Firstly, the sum over all possible decompositions 
constitutes a sum over all possible products of subdiagrams. The kinematic factors may be 
combined into a single term, given that the ^(H) factors commute. Each such product also 
contains a string of exponentiated colour factors, one for each subgraph (i.e. ran = if 
graph H is not present). All permutations of the exponentiated colour factors are present, 
however one should not multiply count permutations obtained by interchanging colour 
factors corresponding to the same subdiagram. Thus, one may write the colour part as a 
sum over all permutations, but with inverse factorial factors m#! to compensate for the 
overcounting. 

Any product of kinematic factors can always be expressed as a sum of kinematic factors 
of graphs D, whose topology is consistent with the diagrams which enter the product. More 
formally, one may write 

HT(Hr« =Y,T(D)N Dl{mH} , (4.11) 

H D 

where -/Vi)|{ mH } is the number of ways in which the diagram D can be obtained from the 
product of kinematic factors in the decomposition given by {ran} (note that this may be 
zero, for the cases in which D cannot be formed from the decomposition). Substituting 
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this into eq. (4.10), the exponentiated amplitude becomes 

ex P fj2m)c(H)\ = f;l e (iw) 

l H ) n=0 ' {m ff } D \ H / 

x p(#i) mi <5(tf 2 ) m2 . . . + perms] . (4.12) 



Currently in the sum over decompositions, it is understood that the sum of all mu is equal 
to n. One may thus replace the multiple sums over n and {m^} with a single sum over all 
decompositions to obtain 

exp(E^)^)l=E^)E^ 1 fn^ ! ) 1 

I H ) D {m H } ' \ H I 

x p(#i) mi C(tf 2 ) m2 . . . + perms] , (4.13) 

where we have also pulled the sum over D to the front. 

Eq. (4.13) must be equal to the original amplitude, which we may express as a sum 
over all diagrams: 

J2?(D)C(D), (4.14) 

D 

where C(D) is the conventional colour factor of graph D. Equating coefficients of T{D) 
with eq. (4.13), one finally finds 

C(D)= E N ° l{ ™ H} (llm H \\ p(^) mi C(F 2 ) m2 ... + perms] , (4.15) 

{m H } n ' \H J 

with n = ^2 H m,H- This formula expresses the conventional colour factors in terms of the 
exponentiated colour factors, and thus is the generalisation of Gatheral's formula to the 
multiparton case. A couple of examples for the application of eq. (4.15) are provided in 
appendix B. The same examples will be analysed in the next section using the replica trick 
method, as well as by direct exponentiation. 

Note that if the colour factors C(D) commute with each other (as happens for two 
eikonal lines), the colour factor of eq. (4.15) becomes 

C , (i^i) mi C(F 2 ) m2 . . . + perms = n\ ]\ C(H) m «, (4.16) 

H 

so that one has 



{rn„} H 

This is precisely the result of Gatheral's formula, albeit expressed in the notation of [89]. 
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5. Exponentiated Colour Factors in special cases 

In the previous sections, we have set up the replica trick for multiparton scattering, and 
presented a general algorithm for calculating the exponentiated colour factor of any graph. 
Furthermore, we have used the algorithm to derive explicit combinatoric formulae for expo- 
nentiated colour factors, which can in principle be applied to arbitrary order in perturbation 
theory. In this section, we present particular cases of applying these results, and our aim 
is twofold. Firstly, we wish to illustrate and thus clarify the somewhat technical discussion 
of the previous sections. Secondly, we will see structures emerging in the exponent which 
allow us to meaningfully generalise the notion of webs from two eikonal line scattering to 
the L-parton case. 

As already stated in the introduction, in the L-parton case one is led to abandon the 
assumption that webs must be single diagrams. Rather, one finds closed sets of diagrams in 
the exponent, whose members are related by permutations of gluons on the eikonal lines, 
and mix only with each other. These diagrams include, in general, ones with reducible 
topologies, with subgraphs that are separable by cutting the eikonal lines, as shown for 
example in figure 9. Such diagrams have sub divergences, in contrast with the familiar 




Figure 9: Example of a reducible topology in which three partons are linked by soft gluon ex- 
changes. Such diagrams have subdivergences, and yet they do contribute to the exponent (as part 
of a multi-diagram web). 

two-eikonal line irreducible webs. 

Before considering this general structure, it is worth noting that there is one way by 
which the irreducibility criterion of webs remains relevant to the multiparton case. This is 
the fact that diagrams that are reducible by mere separation at the hard vertex - namely 
those containing separate clusters of partons that do not communicate by gluon exchanges 
- never contribute to the exponent 11 . We examine this in the following subsection. 

5.1 Diagrams with distinct parton clusters 

As a first demonstration of the replica trick, we consider diagrams in which distinct clusters 
of partons occur, each of which is internally linked by soft gluons, but between which there 

11 Another case in which only irreducible diagrams contribute to the exponent is the planar (large 't Hooft 
coupling) limit, as discussed already, albeit from a different point of view, in [90]. 



- 26 - 



are no exchanges. As an example, we consider explicitly the form shown in figure 10, and 
demonstrate that such diagrams do not contribute to the exponent, independently of the 
internal structure of the subdiagrams G and H. The proof clearly generalises to more than 
two clusters. 




Figure 10: Example of a diagram in which distinct clusters of partons are linked by soft gluon 
exchanges. 



The contribution to the exponentiated colour factor from a given assignment of replica 
numbers in figure 10 is 



N 



(5.1) 



where G{ is the part of subdiagram G composed of replica number i, and C{Gi) its colour 
factor. The matrices C(Gi) and C(Hj) (for any i and j) commute with each other because 
they act on different parton indices. Thus one may rewrite eq. (5.1) as 



A? 



~[{C(G l 



N 



n c(Hj) 



3=1 



(5.2) 



and one sees that the colour factors for subdiagrams G and H are completely independent 
of each other. The number of ways of forming the subdiagram GH is then given by the 
product of the numbers of ways of forming diagrams G and H separately. Given that 
each of the latter is at least O(N), it follows that diagrams of the form shown in figure 10 
are 0(N 2 ), that is their colour factor in the replicated theory, Cn(GH), has no linear 
component in N, so they do not contribute to the exponent. 

5.2 Webs connecting three lines at two loops 

In the previous subsection we have seen an example where reducible diagrams do not con- 
tribute to the exponent. In general, in the multiparton case, things are more complicated 
than this, and in contrast to the two-line case certain reducible diagrams, such as those of 
figure 9, do contribute. We start by examining two-loop diagrams, in which at most three 
partons may be connected. These diagrams have already been studied in the context of 
soft anomalous dimensions [58,68,74], thus we will be able to demonstrate the replica trick 
in a familiar context, and also reinterpret the previously known results in terms of webs. 
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First, there is the diagram in which three gluons are connected by a three-gluon vertex, 
shown in figure 11. It is clear that this diagram is a web from the replica trick, as it is 0{N) 
in the replicated theory (only one replica can be involved due to the lack of self-interactions 
between the replicas). It was indeed shown to be present in the soft anomalous dimension 
matrix for massive eikonal lines [68] (it in fact vanishes for massless eikonal lines, once 
kinematic information is included [58]). 




Figure 11: Example two- loop diagram in which three partons are connected via a three gluon 
vertex. 

There are also the diagrams in which three lines are connected by two single gluon 
exchanges, shown here in figure 12. Using a similar notation to before, we denote by C{a) 
and F{a) the colour and kinematic parts respectively of diagram a (and similarly for b). 
We may then evaluate the exponentiated colour factors using the algorithm defined in 
section 3.2. First one identifies distinct connected pieces, and assigns a replica number to 
each. This has already been carried out in figure 12, where replica numbers i and j are 
introduced. Next, one considers all possible hierarchies h of the replica numbers. There 
are three such hierarchies: i = j , i < j and i > j. For each h, one must evaluate the colour 
factor in the replicated theory (lZ[C(D)\h\) for each diagram D. 




(a) (b) 



Figure 12: Example two-loop diagrams involving two gluon exchanges between three lines. 
Considering i = j, the exponentiated colour contributions from figure 12(a) and (b) 

are 

M N (i = j)K[C(a)\i = j] = NC(a) (5.3) 
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and 

M N (i = j)K[C(b)\i = j] = NC(b). (5.4) 

That is, these are the conventional colour factors for these diagrams, weighted by M^{i = 
j), the number of possible assignments of replica numbers i and j such that i = j, which 
is equal to the number of replicas N. 

Next there is the case i < j. Given the replica ordering in the replicated theory, 
the gluons on the middle line in diagram (a) get interchanged, so that the colour factor 
corresponds to that of diagram (b). However, no such reordering occurs in diagram (b). 
Thus, the contributions from diagrams (a) and (b) are 

M N {t < j)K[C(a)\i < j] = N{N ~ 1] C(b) (5.5) 

and 

M N (i < j)K[C(b)\i < j] = N{N ~ 1] C(b) (5.6) 

respectively, where the prefactor arises from the number of ways of choosing i < j. 

Finally, one has the case i > j. This is similar to the case i < j, except that reordering 
of gluons occurs in diagram (b) but not in (a), and the two diagrams give contributions 

M N (i > j)K[C(a)\i > j] = N{N ~ 1] C(a) (5.7) 

and 

M N (i > j)K[C(b)\i > j] = N{N ~ 1) C(a). (5.8) 

Summing over all hierarchies h (as in eq. (3.15)), the total colour factors of diagrams (a) 
and (b) in the replicated theory are 

N N 2 
C N (a) = - [C(a) - C(b)] + — [C(a) + C(b)] . (5.9) 

N N 2 
C N (b) = - [C(b) - C(a)] + — [C(a) + C(b)] . (5.10) 

The exponentiated colour factors C(a) and C(b) are given by the coefficient of the 0{N) 
parts of these results (eq. (3.16)). Combining with the kinematic factors for each graph, 
the total contribution to the exponent from the diagrams of figure 12 is 

l -[c(a)-C{b))[Ha)-Hb))- (5.11) 

Thus, the only surviving part which exponentiates has the form of a function which is 
antisymmetric in both colour and kinematics. This agrees exactly with the results of [58], 
which state that the sum of diagrams (a) and (b) can be decomposed into a term which 
is symmetric in both colour and kinematics, and a part which is antisymmetric in both. 
The former can be shown to be a product of lower order diagrams, thus does not appear in 
the exponent of the eikonal scattering amplitude. Here we see this explicitly in eqs. (5.9) 
and (5.10), where the symmetric combination is 0{N 2 ). The antisymmetric combination 
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indeed survives (is O(N)) 12 . This is our first example of the fact that reducible diagrams 
can appear in the exponent. Furthermore, the diagrams of figure 12 form a closed set, 
mixing only with each other in colour space. We will see more examples in the following 
section, where we consider diagrams at three loop order. 

5.3 Webs in multi-parton scattering: three loops 

Having shown examples of applying the replica trick to diagrams that have been previously 
studied using other methods, we now turn our attention to cases that have not been 
previously explored, at three-loop order. Although we shall consider here the specific case 
of three loops (and, in the next section, one set of four-loop diagrams) we emphasise that 
the replica trick algorithm, as set up in section 3, can be directly applied at any order in 
the exponent: it does not require prior analysis of lower-order graphs. 

Further to providing examples for using the replica trick algorithm, we will explicitly 
check in each case that the linear combinations of colour factors found from the replica 
trick agree precisely with what is obtained upon calculating the exponent by subtracting 
the exponentiation of lower-order diagrams from the full result. The latter calculation is 
rather long, and it requires manipulating products of colour factors well as kinematic parts 
of lower-order diagrams. This comparison therefore illustrates how effective the replica 
trick algorithm is in determining directly the structure of the exponent. 

Replica trick computation of PF(i 2,2,i) 

We first focus on the set of diagrams shown in figure 13, where we have labelled each 
diagram as shown in the figure. Considering first diagram (3a), one must again assign 
replica numbers to the distinct connected pieces, and consider all possible hierarchies h. 
Labelling the upper, right-hand and lower gluons by i, j and k respectively, the list of 
hierarchies is given in table 1. The second column shows, for each h, the replica ordered 
colour factor (corresponding to the conventional colour factor of one of the graphs in 
figure 13). The third column shows the multiplicity factor M^{h) i.e. the number of ways 
of choosing the number of replica numbers i, j and k out of N, given h. Finally, the fourth 
column shows the O(N) part of this multiplicity factor, which enters the expression for the 
exponentiated colour factor (eq. (3.16)). 

Adding things together, the 0(N) part of the colour factor for this diagram in the 
replicated theory (i.e. the exponentiated colour factor in the unreplicated theory) is 

C(3a) = * [C(3a) - C(36) - C(3c) + C(3d)\ . (5.12) 

One may evaluate the exponentiated colour factors of the remaining diagrams in figure 13 

12 Note that the antisymmetric combination of kinematic factors in eq. (5.11) actually vanishes in the case 
of massless eikonal lines [58], although not in the case of massive lines [68]. 
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[[1],[1,2],[3,2],[3]]=(3a) [[1],[2,1],[3,2],[3]]=(3b) [[1 ],[1 ,2],[2,3],[3]]=(3c) [[1 ],[2,1 ],[2,3],[3]]=(3d) 

Figure 13: Example three loop diagrams. 



h 


K[(3a)\h] 


Mjv(/i) 


©(iV 1 ) part of M N (h) 


i = j = k 


C(3o) 


TV 


1 


i = j,k > i 


C(3a) 


AT(7V - l)/2 


l 

2 


i = j, k < i 


C(3c) 


AT(7V - l)/2 


1 
2 


i = k, j > i 


C(3d) 


AT(7V - l)/2 


1 
2 


i = k, j < i 


C(3o) 


AT(iV - l)/2 


1 
2 


j = k,i> j 


C(3o) 


AT(jV - l)/2 


1 
2 


j = k,i < j 


C(36) 


AT(7V - l)/2 


1 
2 


i < j < k 


C(36) 


7V(iV- l)(JV-2)/6 


1 

3 


i < k < j 


C(3d) 


iV(i\T — 1)(JV — 2)/6 


1 

3 


j < i < k 


C(3a) 


iV(iV- l)(JV-2)/6 


1 

3 


j < k < i 


C(3o) 


N(N- l)(JV-2)/6 


1 

3 


k < i < j 


C(3d) 


7V(JV- l)(iV- 2)/6 


1 

3 


k < j < i 


C(3c) 


iV(JV — l)(JV-2)/6 


1 

3 



Table 1: Replica trick analysis for diagram (3a), as shown in figure 13. The replica indices i, j 
and k label the gluons on the upper, right-hand and lower parts of the diagram respectively. 

in a similar fashion to obtain 

C(3b) = \ \-C{3a) + C(3b) + C(3c) - C(3d)} (5.13) 

C(3c) = - [-C(3a) + C(3b) + C(3c) - C(3d)] (5.14) 

C(3d) = - \C(3a) - C(3b) - C(3c) + C(3d)} , (5.15) 
6 

so that the total contribution to the exponent from these diagrams is 

W(i i2 ,2,i) = g(c(3a)-C(36)-C(3c) + C(3d)) (r(3a)- 2T(3b) -2J r (3c)+J r (3d)) , (5.16) 

where, as usual, J~(D) denotes the kinematic part of diagram D, and on the left-hand side 
we have used the notation of eq. (1.4) denoting the number of gluon attachments on each 
parton line. 
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It is also instructive to see the above results derived from the closed form solution for 
exponentiated colour factors of eq. (4.3). We carry out this calculation in Appendix A. 
Similarly to the two-loop case considered in the previous section, we see that the diagrams 
of figure 13 form a "closed set" of diagrams which mix with each other in the exponent. 
Equation (5.16) is thus the result for the "web" formed from these diagrams. Furthermore, 
a non-trivial degree of combinatorics enters eq. (5.16). One might have not anticipated, 
for example, the factors of two in this expression. We shall return to this point below. 

Replica trick computation of W^ 2 ^,i) 

First let us consider another example of a closed set of diagrams, namely the set (2, 3, 1) 
in the notation of eq. (1.4). These are shown in figure 14, where we also introduce some 
convenient labels. 






[[1,2],[2,1,3],[3]]=(3A) 



[[1,2], [2,3,1], [3]]=(3B) 



[[1,2], [3,2,1], [3]]=(3C) 






[[1 ,2],[1 ,2,3],[3]]=(3D) [[1 ,2],[1 ,3,2],[3]]=(3E) [[1 ,2], [3,1 ,2],[3]]=(3F) 

Figure 14: Example diagrams in which three parton lines are linked by three gluon exchanges. 



Using the replica trick as outlined above, one may show that the exponentiated colour 
factors are 

1 



C(3A) 
C(3B) 
C{3C) 
C(3D) 
C(3E) 
C(3F) 



3C{3A) - 2C(3D) - 2C(3E) - 3C(3C) + 4C(3F) 

- 3C(3A) + 6C{3B) + C(3D) - 2C{3E) + C(3F) - 3C(3C) 
3C(3C) - 2C{3E) - 2C(3F) - 3C{3A) + AC{3D) 
C(3D) - 2C(3E) + C{3F) 

- C{3D) + 2C{3E) - C(3F) 
C(3D) - 2C(3E) + C{3F) 



(5.17) 
(5.18) 
(5.19) 
(5.20) 
(5.21) 
(5.22) 
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such that the total contribution to the exponent from these diagrams is 



W, 



(2,3,1) 



+ 



1 

6 

l r 



+ 



+ 



l r 

6 
1 



+ 



3 L 



3C(3A) - 3C(3C) - 2C(3D) - 2C(3E) + AC{3F) F{3A) 

- 3C(3A) + 6C{3B) - 3C{3C) + C{3D) - 2C{3E) + C(3F)]f(3B) 

- 3C(3A) + 3C(3C) + 4C(3D) - 2C(3E) - 2C(3F)1 F(3C) 

C{3D) - 2C(3E) + C{3F)\f(3D) 

C{3D) + 2C{3E) - C{3F) F{3E) + ^ C(3D) - 2C(3E) + C(3F) F(3F). 

(5.23) 



Once again, it is instructive to rederive the above results using eq. (4.3) - see appendix A. 
Here the combinatorics of how each coefficient is obtained are even more complicated, and it 
is useful to reproduce the above results by explicitly exponentiating lower-order diagrams. 



Calculation of W(i 5 2,2,i) by direct exponentiation of lower orders 

Consider again the diagrams of figure 13, whose contribution to the exponent, as found 
from the replica trick, is given by eq. (5.16). One may verify this result explicitly as 
follows. One exponentiates all one and two loop diagrams (with the appropriate exponen- 
tiated colour factors) that lead to the diagrams whose form matches those of figure 13. 
Then, one subtracts the result from the appropriate contribution to the unexponentiated 
amplitude. The result is, by definition, the contribution to the exponent from these three 
loop diagrams. 

The aim of this exercise is both to demonstrate how exponentiation of lower-order 
diagrams proceeds in the multileg case, and to highlight the drastic simplification achieved 
by the replica trick. We will see that the results we have obtained above, even though they 
relate to three-loop diagrams of relatively simple topology, are much more complicated 
to derive using explicit exponentiation, which involves manipulating both kinematic and 
colour factors of all relevant lower-order diagrams. The reader who is not interested in 
seeing such a rederivation of the above results may skip to the end of this section, where 
we discuss the general structure of our results in more detail. 

To be more specific, one may write the eikonal scattering amplitude in either conven- 
tional or exponentiated form as 



A = Y,A n a n s = exp(^5 m cc) , 

n \ m / 



(5.24) 



where the coefficients A, and Bi are matrix-valued in the space of possible colour flows, 
thus do not commute in general. In the notation introduced above, the {Ai} are given by 



A n = Y J C{G)F{G), 

G 



(5.25) 
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where the sum is over all Feynman diagrams G contributing at C(a"). By expanding the 
exponential on the right-hand side, the {Bi} and {Ai} can be related to each other. Up to 
three-loop order, one has 

A 1 = B 1 

A 2 = B 2 + X -B\ 

A 3 = B 3 + ^{B 1 ,B 2 } + ^Bf. (5.26) 

In order to rederive the replica trick result (i.e. that B% is given by eq. (5.16)) using (5.26), 
one thus needs to evaluate the one and two-loop coefficients B\ and B 2 , and subtract 
the appropriate combination from A3. In calculating B\ and B 2 , one need only consider 
diagrams that, when exponentiated, will give rise to three loop diagrams of the form shown 
in figure 13. The relevant diagrams are shown in figures 15 and 16. We have seen some of 
these diagrams already. However, we redraw them all together here for convenience, and 
also to introduce labels that will be used in what follows. From previous results we know 
that 

Bx = C(la)F(la) + C(lb)F(lb) + C(lc)F(lc) 

B 2 = \{c(2a) - C(2b)) (>(2a) - F(2b)) + \(c(2c) - C(2d)) (f(2c) - T(2d)), (5.27) 

where we have kept only those terms which contribute to eq. (5.16). The first result follows 
immediately from B± = A\ , and the second result was first derived in [58] (we have also seen 
part of this result in eq. (5.11)). In combining these results to explicitly derive S3, one needs 
to know how to deal with products of eikonal diagrams (e.g. C(la)C(2a)F(la)F(2a)). To 
illustrate the method, it is simpler to first rederive B 2 by exponentiating the one-loop 
result. 




[[i].[i].[].[]]=(ia) m,[i],[i],[]]=(ib) [[],[], [1],[1]]=(1c) 

Figure 15: One loop diagrams which, when exponentiated, contribute to eq. (5.16). 



At two loops (5.26) yields 



B 2 = A 2 - ^Bj, (5.28) 



where 



A 2 = C(2a)F(2a) + C (2b) J 7 (2b) + C(2c)F(2c) + C(2d)F(2d) + CWJ^e). (5.29) 



[[1],[1,2],[2],[]]=(2a) [[1],[2,1],[2],[]]=(2b) [[],[1],[1,2],[2]]=(2c) [[],[1 ],[2, 1 ],[2]]=(2d) 

Figure 16: Two loop diagrams which, when exponentiated, contribute to eq. (5.16). 




[[1],[1],[2],[2]]=(2e) 

Figure 17: Two- loop diagram present in the amplitude, but absent from the exponent. 



The final term contains the diagram shown in figure 17 which, as we will see, does not 
appear in the exponent. 

Evaluating B\ from eq. (5.27) gives 



l R2 _lr 



7(la)7(lb) {C(la),C(lb)} + JT(la)jr(l c ) {C(la), C(lc)} 
+ T(lb)T(lc) {C(lb), C(lc)} + . . . 1 . (5.30) 



On the right-hand side, we have kept only products of diagrams which give rise to graphs 
of the same form as those shown in figure 16. Furthermore, we have used the fact that the 
kinematic factors commute ([J 7 (G) , J 7 (H)] = 0), but the colour factors do not in general. 
Considering the first term, the product of kinematic factors gives 

/'OO /'OO /'OO poo 

J 7 (la) J 7 (lb) = V12V23 / dh / dt 2 / ds 2 / ds 3 D(tiVi - t 2 v 2 )D(s 2 v 2 - s 3 v 3 ), 

JO JO JO JO 

(5.31) 

where vi is the 4- velocity (or rescaled 4-momentum) of the i th eikonal line, and where we 
introduce the notation v ij = Vi ■ Vj for the scalar product of Vi and Vj. Here D(x — y) is the 
position-space propagator for a soft gluon propagating between x and y. We have ignored 
coupling factors etc. for brevity. The colour factors may be given in Catani-Seymour 
notation [114] as 

C(la) = T?T£, C(lb) = T^T|, (5.32) 
(with T" a colour generator in the appropriate representation of parton i) so that 

{C(la),C(lb)} = TUT a 2 ,T b 2 }T b 3 . (5.33) 
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Then the first term on the right-hand side of eq. (5.30) gives 

; poo poo poo poo 

^v 12 v 23 T a l {T^,T b }T b dh dt 2 ds 2 ds 3 D(t lVl - t 2 v 2 )D(s 2 v 2 
1 Jo Jo Jo Jo 



S3V 3 ). 

(5.34) 



Similarly, the second and third terms give 

/'CO f'OO /'CO' /'CO' 

v i2 v 23 T^T c 3 T c 4 dh dt 2 dt 3 dt A D{t lVl -t 2 v 2 )D(t 3 v 3 -Uv 4 ) (5.35) 
Jo Jo Jo Jo 



and 



j />00 /"OO /'CXD rOO 

-v 12 v 23 T b {T b 3 ,T c 3 }T c 4 ds 2 ds 3 dt 3 dt 4 D{s 2 v 2 -s 3 v 3 )D(t 3 v 3 -t 4 v 4 ) 
1 Jo Jo Jo Jo 

(5.36) 

respectively. From eq. (5.29), one has 

/'OO poo pt 2 poo 

A 2 =v 12 v 23 T^T&T* / dh dt 2 ds 2 ds 3 D(t lVl - t 2 v 2 )D(s 2 v 2 - s 3 v 3 ) 

Jo Jo Jo Jo 

poo poo poo poo 

TfT^Tl / dh / dt 2 / ds 2 / ds 3 D{t lVl - t 2 v 2 )D(s 2 v 2 - s 3 v 3 ) 
Jo Jo Jt 2 Jo 

poo poo poo poo 

T b 2 T b 3 T c 3 Tl / ds 2 / dt 3 / ds 3 / dt A D(s 2 v 2 - s 3 v 3 )D(t 3 v 3 - t 4 v 4 ) 
Jo Jo Jt 3 Jo 

poo poo ptz poo 

+ T^T|T2 / ds 2 / dt 3 / ds 3 / dt 4 D(s 2 v 2 - s 3 v 3 )D(t 3 v 3 - t 4 v 4 ) 
Jo Jo Jo Jo 

poo poo poo poo 

+ T?Tf PgT^ / dh / dt 2 / dt 3 / dt 4 D(hv! - t 2 v 2 )D{t 3 v 3 - t 4 v 4 ) , (5.37) 
Jo Jo Jo Jo 

where the five terms correspond to diagrams (2a)-(2e) respectively. Using eqs. (5.28,5.34- 
5.36) one finds 



+ ' 
+ ' 



poo poo poo 

B 2 = v 12 v 23 \-Tt\T^T b 2 \T b 3 dh J dt 2 J ds 3 



t 2 poo 

ds 2 — I ds 2 
o Jt 2 

oo />oo 



poo poo poo 

D(hvi - t 2 v 2 )D(s 2 v 2 - s 3 v 3 ) + -T b [T b , T C 3 ]T C 4 / ds 2 / dt 3 / dt 4 

1 Jo Jo Jo 

f'OO rtz 

/ ds 3 - c 
Ju Jo 



D(s 2 v 2 - s 3 v 3 )D(t 3 v 3 - t 4 v 4 ) 

Ut 3 JO J 

= l - (C(2o) - C(26)) (.F(2a) - F(2b)) + 1 (C(2c) - C(2d)) (.F(2c) - .F(2d)) , (5.38) 

which indeed agrees with eq. (5.27). 

Having seen explicitly how exponentiation of a subset of diagrams works at two loops, 
we now consider the three loop example of figure 13 by evaluating B 3 , given by eq. (5.26) 
as 

i i 

(5.39) 



B 3 = A 3 - l -B\- l -{B ll B 2 } 1 



where 



A 3 = C(3a)T(3a) + C (3b) J 7 (3b) + C(3c)F(3c) + C{3d)F(3d). 



(5.40) 



- 36 - 



Firstly, one has 
\ B l = I (c{la)F{la) + C{lb)F{lb) + C{lc)F{lc)f 

= * F(la)F{lb)F{lc) (c(la)C(lb)C(lc) + C{la)C{\c)C(lb) + C(lb)C(la)C(lc) 
+ C(lb)C(lc)C(la) + C(lc)C(lo)C(16) + C(lc)C(16)C(la)) + • • • , (5.41) 

where we have kept only combinations in the cubic product that lead to diagrams of the 
form shown in figure 13, and have also used the fact that the kinematic factors commute. 
Considering the first term in the sum over products of colour factors, this is 

C(la)C(lb)C(lc) = (T?T2)(T§T|)(T§T2) = T^T^T^ 

= C(3c). (5.42) 

Note we have used the rule in combining the colour factors that in each product of matrices 
associated with the same parton line (e.g. T^Tg), the matrices are ordered from the external 
line towards the hard interaction. Evaluating the remaining terms in eq. (5.41), one finds 

= 1 j-(i a ) j-(ife) j-(i c ) (c(3c) + 2C(3a) + 2C(3d) + C(36)) . (5.43) 

The kinematic prefactor has the explicit form 

roc poo poo poo poo 

F(la)F(lb)F(lc) = vi2?>23f34 / dti / dt2 / dt^ I ds2 / ds^ 

Jo Jo Jo Jo Jo 

poo 

x / dt 4 D(tiVi - t 2 v 2 )D(s2V2 - s 3 v 3 )D(t 3 v 3 - t 4 v 4 ). (5.44) 
Jo 



All of the kinematic factors from now on will involve the same integrals over {U}, and the 
same product of propagators and Vij factors. Thus we may shorten notation by writing 
each kinematic factor according to its s 2 and S3 integrals i.e. 

poo poo 

F(la)F(lb)F(lc) = / ds 2 / ds 3 . (5.45) 
Jo Jo 

We may express this in terms of the kinematic factors of the individual diagrams in figure 13. 
After writing 

00 pti poo 

dsi = dsi+ I dsi, (5.46) 
Jo Ju 

eq. (5.45) becomes 

F{la)F{lb)F{lc) = F(3a) + F{3b) + F(3c) + F(3d), (5.47) 
where we have identified 

rt 2 pts poo pUi 

'*3 



F(3a) = / ds2 I ds 3 , F(3b) = / ds2 / ds 3 
Jo Jo Jt 2 Jo 

pt'2 poo poo poo 

F(3c) = ds 2 ds 3 , F{3d) = / ds 2 / ds 3 . (5.48) 
Jo Jt A Jt 2 Jt z 
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Thus one has 

^Bf = ^ (>(3a) + F(3b) + JF(3c) + JF(3d)) (2C(3a) + C(36) + C(3c) + 2C*(3d)) . (5.49) 

This result is already interesting, in that some nontrivial combinatoric factors appear. The 
three-loop diagrams do not occur equally, with a separation between diagrams (3a, 3d) and 
(3b, 3c). We will see that the same separation persists in the final result for B 3 (as has 
already been found using the replica trick). 
The last ingredient in (5.39) is the term 

^{B U B 2 } = |-F(la)(-F(2c) - T(2dj) \c(la)(c(2c) - C(2d)) + (c(2c) - C(2df)C(la) 

+ ^(lc) (V(2a) - J-(2fe)) \c(lc) (c(2a) - C(26)) + (c(2a) - C(26)) C(lc) 

+ ... (5.50) 

Again we have only kept terms leading to diagrams having the forms shown in figure 13. 
The kinematic factors may be simplified as follows: 

F(la) (T(2c) - F(2d)) = ds 2 ds 3 - ds 2 ds 3 



I* ds 2 [ 

JO Jt: : 

f'OO f' 

- ds 2 
hi Jo 



t :i JO JO 

00 /*oo roo 

dss + / ds2 dsz 

t 3 

ds 3 



rt2 

/ ds 2 / ds 3 
JO JO 



= F(3c) + F(3d) - F(3a) - F(3b); 



(5.51) 



F(lc) (F(2a) - F(2b)) 



t2 />00 /'OO />oo 

ds 2 I ds 3 — I ds 2 I ds 3 
Jo Jt 2 Jo 

ds 3 + I ds 2 I ds 3 

t 3 Jt 2 



t>t2 rt 3 rt2 roo poo rtz 

/ ds 2 / ds 3 + I ds 2 / ds 3 - ds 2 / ds 3 
Jo Jo Jo Jt s Jt 2 Jo 



ds 2 / ds 3 

t 2 Jt z 



= F(3a) + F(3c) - F(3b) - T(3d) . (5.52) 
Furthermore, the colour factors can be evaluated to give 

C(la)(c(2c)-C(2d)j + (c(2c)-C(2d)^C(la) = C (3c) - C (3a) + C (3d) - C (3b) . (5.53) 
and 

C(lc)(c(2a)-C(2b)^j + (c(2a)-C(2b)^C(lc) = C '(3c) - C (3d) + C '(3a) - C (3b) . (5.54) 
Putting things together, one finds 

B 2 } = X - \c(3a) (T(3a) - JF(3d)) + C(3b) (?(3b) - JF(3c)) + C(3c) (>(3c) - F(3b)^j 



+ C(3d)[F(3d)-F(3a) 



(5.55) 
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Finally, combining eqs. (5.39, 5.40, 5.49, 5.55) one finds 




which agrees precisely with the replica trick result of eq. (5.16). Furthermore, we have 
learnt the following lessons: 

• The exponentiation of lower-order diagrams results in products of colour factors, as 
well as products of kinematic factors, which can be recast as linear combinations 
of the colour and kinematic factors of recognisable diagrams. These manipulations 
bring about non-trivial combinatorics. 

• The replica trick greatly simplifies the calculation. The degree of simplification be- 
comes yet more apparent at higher orders. 

Calculation of W(2,s,i) by direct exponentiation of lower orders 

As a further demonstration, one may consider the diagrams of figure 14. The combination 
of these diagrams which appears in the exponent is given in eq. (5.23), and one may check 
that the same result is obtained on subtracting the relevant exponentiated one and two 
loop graphs from the three loop amplitude. In this case, the terms Bi in the exponent up 
to three loops are given by eq. (5.26) with 



where the second line contains the diagrams shown in figure 18. As before, we first calculate 



Figure 18: Two loop diagrams which, when exponentiated, contribute to the diagrams shown in 
figure 14. 

the two loop contribution to the exponent, before extending the calculation to three loops. 



A 3 



C(la)F{la) + C(lb)F(lb) 

C{2a)F{2a) + C(2b)F(2b) + C(2f)F(2f) + C(2g)F(2g) 
C(3A)F(3A) + C(3B) F(3B) + C{3C)F{3C) + C{3D)F(3D) 
+ C(3E)F(3E) + C{3F)F{3F), 



(5.57) 




[[1,2],[1,2],[],[]]=(2f) 



[[1,2],[2,1],[],[]] = (2g) 
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For this set of diagrams, one has 



\b\ = \ (c(la)T(\a) + C(lb)T(lb)) 

= \ C\\a)T\\a) + \{C(la\ C(lb)}F(la)F(lb) + ..., 



(5.58) 



where the ellipsis denotes terms which do not contribute to diagrams of the form shown in 
figure 14. The kinematic factor in the first term is 

pod pod pod pod 

F 2 (la)=v 2 2 I dti / dt 2 / dsi / ds 2 D(t 1 vi - t 2 v 2 )D(s\Vi - s 2 v 2 ), (5.59) 
Jo Jo Jo Jo 

which must be rewritten in terms of the kinematic parts of the diagrams shown in figure 18. 
We parametrize the integrals along v\ by t\ and s±, the ones along v 2 by t 2 , s 2 and S3, and 
the one along v 3 by £3. Given these notations, we note that the integrals considered share 
the same ti integrals as well as the same propagator factors, while the limits on the s±, s 2 
and S3 integrations differ between them. We can therefore represent the kinematic part of 
each diagram by its Sj integrals 13 i.e. for the two-loop integrals in figure 18, 



F{2f) = dsi ds 2 or / dsi / 

Jtx Jt 2 Jo Jo 

J 7 (2g) = ds\ ds 2 or / dsi \ 

Jti Jo Jo Jt 2 



ds 2 



ds 2 



(5.60) 



Note that each diagram can be rewritten in two ways, which correspond to different la- 
bellings of gluons in terms of the parameters Sj and ij. One then has 



F 2 {la) = 



rti roc 1 r rt 2 r 

/ dsi+ dsi / ds 2 + I 

Jo Jt x J YJo Jt : 

2(j-(2/)+J-(2 5 )). 



(5.61) 



One also needs the kinematic combination 

POO POD POD 



POD POD POD POD 

J 7 (la) J 7 (lb) = v 12 v 2 3 / dh / dt 2 I ds 3 / dhD(tiVi - t 2 v 2 ) D(s 3 v 2 - t 3 v 3 ) 
Jo Jo Jo Jo 

dt 3 



OD roo 

dU 



vi 2 v 23 I dti 

10 Jo 
»t 2 roo 
ds 3 + / ds 3 
lt 2 



D(tivi - t 2 v 2 ) D(s 3 v 2 - t 3 v 3 ) 



[ ds 3 + f 
JO Jt : 

= F(2a) + F(2b). 

The colour factors appearing in eq. (5.58) combine to give 

C 2 (la) = C(2f) 
{C(la), C(lb)} = C(2a) + C(2b). 



(5.62) 



(5.63) 



13 Similar notation will be used below at three loops. 
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Putting things together, one finds 



^ = C(2f) (F(2f) + F{2g)) + \ (c(2a) + C(26)) (>(2a) + ^(26)) . (5.64) 

Combining this with eqs. (5.26, 5.57), the two-loop contribution to the exponent for the 
diagrams of figure 14 is 

B 2 = (C(2g) - C{2f))F(2g) + \ (c(2a) - C{2b)) (f{2o) - F{2b)) , (5.65) 

which agrees with the results of the replica trick, and also those of [58] . 

Next, we calculate the three loop contribution to the exponent from the diagrams of 
figure 14. From eq. (5.26), one first has the cube of the one loop exponent 

^= l -(c(la)F(la) + C{lb)F{lb)f 

= - F 2 (la) F(lb) (c 2 (la)C(lb) + C(la)C(lb)C(la) + C (lb)C 2 {la^j + ..., (5.66) 

where we again neglect terms corresponding to diagrams which do not have the form of 
figure 14. The kinematic factor must be rewritten in terms of the kinematic factors of each 
diagram. The kinematic factor for diagram (3A) is 

poo poo poo pti poo pt2 

F(3A) = v 2 2 v 23 / dti / dt 2 / dt 3 / dsi / ds 2 / ds 3 
Jo Jo Jo Jo Jt 2 Jo 

x D(t\vi - t2V 2 )D(siVi - s 2 v 2 )D(s3V 2 - t 3 v 3 .) (5.67) 



Again we may neglect the propagator factors and U integrals, and represent each diagram 
by its Si integrals. The results are 



rti poo pt2 poo pt2 pS2 

^"(3^4) = / ds\ / ds 2 / dss or / ds\ I ds 2 I ds^; 
Jo Jt 2 Jo Jti Jo Jo 

pti poo pS2 poo pt 2 pt2 

F(3B) = / ds\ / ds 2 / ds3 or / ds\ / ds 2 \ ds^; 
Jo Jt 2 Jt 2 Jti Jo Js 2 

pt\ poo poo poo pt2 poo 

F(3C) = / dsi / ds 2 / ds3 or / dsi / ds 2 / ds% 

Jo Jt 2 Js 2 Jti Jo Jt 2 

pti pt2 pS2 poo poo pt2 

F(3D) = / ds\ / ds 2 / ds3 or / ds\ / ds 2 / ds3; 
Jo Jo Jo Jti Jt 2 Jo 

pti pt2 pt2 poo poo pS2 

F(3E) = / ds\ / ds 2 / ds3 or / ds\ / ds 2 / ds3 

Jo Jo J s 2 Jti Jt 2 Jt 2 

pti pt2 poo poo poo poo 

F(3F) = dsi ds 2 / ds 3 or / ds 1 / ds 2 / ds 3 . (5.68) 

Jo Jo Jt 2 Jti Jt 2 J S2 
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where we have used the fact that each diagram may be written in different, but equivalent 
ways. Then one can write 

^oo roo roo 

J- 2 (la)J-(16) = / dsi / ds 2 / ds 3 
Jo Jo Jo 

rti rt 2 r rS2 rti roo n roo rt 2 r />s 2 

= / dsi ds2 I ds 3 + J ds 3 + / ds 3 + dsi / ds 2 / ^3 
jo jo L.J0 Js 2 Jt 2 J Jti jo Uo 

rt\ roo r /><2 /<S2 /-oo 

/ dsi / ds2 I ds 3 + I ds 3 + / ds 3 
Jo Jt 2 Uo Jt 2 Js 2 



(5.69) 



rt 2 roo 

+ ds 3 + dt 

J s 2 Jt 2 

roo roo r rt 2 rs 2 roo 

+ / dsi ds 2 / ds 3 + ds 3 + / (is 
Jti Jt 2 Uo Jt 2 Js 2 

= 2(T(3A) + T(3B) + JT(3C) + F(3D) + JF(3£) + T(3F)). 
The colour factors in eq. (5.66) give 

C 2 (la)C{lb) = C{3D); 
C(lo)C(16)C(lo) = C(3£); 
C(16)C 2 (la) = C{3F), 

so that one has 

^ = J (C(3D) + C{3E) + C{3F)){T{3A) + F{3B) + T(3C) + F{3D) + T(3E) + T(3F)). 

(5.71) 

Next, one has the contribution of one and two loop webs from the exponentiation i.e. the 
third term in eq. (5.39), which in the present case gives 

* B 2 } = \{C{la), C{2a) - C{2b)}F(la)(F(2a) - F(2b)) 



(5.70) 



+ l -{C{lb),C(2g) - C{2f)}F{\b)F{2g). 
The first kinematic factor on the right-hand side is given by 



(5.72) 



F(la)(T(2a) - T{2b)) = I ds ± I ds 2 I ds 3 - I dsi I ds 2 I ds 3 

'0 JO JO JO JO Js 2 

rt\ rt 2 roo rt 2 rt\ roo 

I ds\ I ds 2 + / ds\ I ds 2 + / ds\ I ds 2 
Jo Jo Jti Jo Jo Jt 2 

f 

J So 



roo roo 

+ / ds\ / ds 

Jt\ Jt 2 



ds 3 - I ds 3 

s 2 

2{F{3A) - F{3C) + F{3D) - F{3F)). 



(5.73) 



Also, one has 



pt\ fOG f'OO 

F{lb)T(2g) = ds! ds 2 / ds 3 
Jo Jt 2 Jo 

rt\ roo r rt 2 rs 2 roo 

= dsi ds 2 I ds 3 + I ds 3 + / ds 3 

Jo Jt 2 Uo Jt 2 Js 2 

= F{3A) + F{3B) + F{3C). 



(5.74) 
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The colour factors in eq. (5.72) give 



{C(la), C(2a) - C(2b)} = C{3D) - C(3F) 

{C(lb), C{2g) - C(2/)} = C(3A) + C(3C) - C(3D) - C(3F). 



(5.75) 



Putting things together, one has 

l -{B l ,B 2 ] = l - (C(3D) - C(3F)) (t(3A) - F{3C) + F{3D) - F{3F)) 

+ X - (C{3A) + C{3C) - C(3D) - C(3F)) (V(3A) + F(3B) + J"(3C)) . 

(5.76) 

Finally, combining eqs. (5.39, 5.57, 5.71, 5.76), one finds that B% is given by exactly the 
same expression as eq. (5.23), thus agreeing with what has already been found using the 
replica trick. 

Discussion 

We have now seen for two non-trivial examples that the replica trick produces the same 
results as can be obtained from the explicit exponentiation of lower order diagrams, as it 
must do. We see that in general, the exponent contains distinct sets of diagrams, related by 
permutations of the gluons on each eikonal line. That such diagrams mix is an expression 
of the fact that they are colour linked. We may then introduce a matrix notation, and 
write the results of eqs. (5.16) and (5.23) as 



(1,2,2,1) 



( ^( 3fl ) ^ 


T 


/ 1 


-1 


-1 


l\ 


/ C(3a) \ 


F{3b) 


1 


-2 


2 


2 


-2 


C{3b) 


F{3c) 


6 


-2 


2 


2 


-2 


C(3c) 


\F{3d)) 




V 1 


-1 


-1 


V 


\C(3d)) 



(5.77) 



and 



W, 



(2,3,1) 



/ F{3A) \ 


T 


/ 


3 


-3 


-2 


-2 


4\ 


/ C(3A) \ 


F{3B) 






-3 6 


-3 


1 


-2 


1 


C(3B) 


F{3C) 


1 




-3 


3 


4 


-2 


-2 


C(3C) 


F{3D) 


6 










1 


-2 


1 


C{3D) 


F{3E) 












-2 


4 


-2 


C(3E) 






V 








1 


-2 


V 


\C(3F)J 



(5.78) 



respectively. In each case a square matrix appears which compactly encodes the mixing of 
each set of diagrams. The general form of such a set is thus 



D D,D> 



(5.79) 



where D and D' label diagrams in a given set, which are mutually related by permutation 
of the gluon attachments to any of the eikonal lines. The set itself is characterized by the 
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number of attachments to each leg (rei, ri2, . . . , n^), as well as by the way subsets of these 
gluons are connected away from the eikonal lines. As usual, C(D) and F{D) denote colour 
and kinematic parts. 

The mixing matrix Rdd' encapsulates the explicit results for the exponentiated colour 
factors in terms of the conventional colour factors. Note that the form of eq. (5.79) also 
applies in the trivial case in which a set contains only one diagram. In such a case, the 
mixing matrix is 1 x 1. The study of web structure then amounts to investigating the 
properties of the mixing matrices Rdd 1 - Having introduced these concepts, we consider 
more examples in appendix C. In the next section, we study general properties of the 
mixing matrices in more detail, starting from a four loop example. 

5.4 Webs in multi-parton scattering: four loops and general discussion 

In this section, we consider the diagrams shown in figure 19, consisting of four parton lines 
connected by four soft gluons. This set of 16 diagrams form a closed set i.e. mix only with 
each other in the exponent. We label these diagrams using the list notation introduced in 
section 3. The mixing Rdd' C(D') for these diagrams is given by 
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[1,2] 


[3,2] 


[4,3], [1,4]] 
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4 


-4 


-2 


-2 


-2 
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[1,2] 


[1,3] 


[4, 3], [4, 2]] 
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-2 
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[1,2] 


[3,2] 


[3,4], [1,4]] 
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2 
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-4 
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[1,2] 


[1,3] 


[3, 4], [4, 2]] 
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-6 
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-4 
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-4 
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[1,2] 


[3,1] 


[4, 3], [4, 2]] 
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-6 
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-4 
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-2 
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[1,2] 


[1,3] 


[4, 3], [2, 4]] 
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-2 
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[1,2] 


[1,3] 


[3,4], [2,4]] 
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[1,2] 


[2,3] 


[4, 3], [1,4]] 
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[1,2] 


[3,1] 


[3, 4], [4, 2]] 
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-2 
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-2 
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-2 
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C| 


[1,2] 


[3,2] 


[4, 3], [4, 1]] 
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-6 


-6 


-6 


-18 


12 


12 


-6 


-18 
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12 


12 


12 


12 
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[1,2] 


[3,1] 


[4, 3], [2, 4]] 
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-6 


-18 


-18 


-18 


-6 


12 


12 


-18 


-6 


-6 


12 


12 


12 


12 





24 ) 




u 


[1,2] 


[2,3] 


[3,4], [4, 1]]; 



(5.80) 



where we have included the vector of colour factors C(D') in order to define the ordering 
of the rows of the matrix. 
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[[1,2], [3,1], [3,4], [2,4]] 



[[1,2], [2,3], [4,3], [4,1]] 



[[1,2], [3,2], [3,4], [4,1]] 



[[1,2], [2,3], [3,4], [1,4]] 





[[1,2], [3,2], [4,3], [1,4]] 



2QQQQQQQQQQQC 
[[1,2],[1,3],[4,3],[4,2]] 





[[1,2], [3,2], [3,4], [1,4]] 



[[1,2], [1,3], [3,4], [4,2]] 







[[1,2], [3,1], [4,3], [4,2]] 



[[1,2], [1,3], [4,3], [2,4]] 



[[1,2], [1,3], [3,4], [2,4]] 



[[1,2], [2,3], [4,3], [1,4]] 





[[1,2], [3,1], [3,4], [4,2]] 



2QQQQQQQQQQQC 

[[1,2],[3,2],[4,3],[4,1]] 





[[1,2], [3,1], [4,3], [2,4]] 



[[1,2], [2,3], [3,4], [4,1]] 



Figure 19: The four loop diagrams which give rise to eq. (5.80). 



As in the three loop cases, there is a non-trivial structure to this matrix. We note 
here some interesting properties. There are two columns of the mixing matrix of eq. (5.80) 
which contain only a single non-zero entry: these are the last two columns. The meaning 
of this becomes clear on examining the last two diagrams in figure 19. We refer to these 
as staircase diagrams, due to the visual similarity with the well-known Escher staircase. 
The last two columns of the matrix then tell us that the colour factors of the staircase 
diagrams do not enter the exponentiated colour factors of the non-staircase diagrams. It 
is easy to understand how this property arises from the replica trick. The replica ordering 
operator TZ can either leave gluon structures intact (in cases where all gluons have the 
same replica number, or are already replica ordered), or disentangle crossed gluons. It 
can never cross gluons that were previously uncrossed, as this results in an inconsistency 
between the replica orderings on different lines. In the present case, one may verify that 
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the staircase diagrams are reordered to non-staircase diagrams if more than one replica 
is present. Furthermore, it is impossible to reorder a non-staircase diagram to make a 
staircase diagram. 

The special feature of those staircase diagrams, which do not enter the exponentiated 
colour factors of other diagrams in the set, is that they have no ultraviolet sub divergences. 
To see this note that in these diagrams the four gluons are maximally linked: one cannot 
shrink any one gluon into the centre of the diagram (the hard vertex) without also pulling in 
all other gluons. This means that these diagrams have a single overall ultraviolet divergence, 
with no subdivergences 14 . It is easy to verify from figure 19 that this is not the case for any 
non-staircase diagram; these all have subdivergences. It follows that the staircase diagrams 
cannot be produced by the exponentiation of lower order diagrams. The ECF of any graph 
consists of the normal colour factor, plus terms which effectively subtract off those parts 
which arise from the exponentiation of lower order graphs. Thus, the colour factors for 
the staircase diagrams cannot enter the exponentiated colour factors of the non-staircase 
diagrams. 

Note that there is also a column containing a single non-zero entry (unit entry) in 
the mixing matrix of eq. (5.78). In this case the diagram corresponding to this column is 
the graph (3B) of figure 14, and one again sees that this is maximally irreducible, in the 
sense that it has no ultraviolet subdivergences (it is not possible to shrink any gluon to 
the centre of the diagram independently of the others). Further examples of diagrams with 
no ultraviolet subdivergences are provided in Appendix C. These are the third diagram 
[[1],[2,1,2],[2],[]] in figure 23, the first and the last diagrams in each of the figures 24 and 26 
and all the diagrams in figure 28 but the second one. In each case there is a corresponding 
column containing a single non-zero entry in the mixing matrix. 

Another remarkable property of the mixing matrices we computed is that they all 
have zero-sum rows, that is the entries in any row of these matrices sum to zero 15 . In the 
notation of eq. (5.79) one has 

J2Rdd>=0, VD. (5.81) 

D' 

This implies a symmetry of eq. (5.79), 

C(D) = Rdd> C(D') = ^ Rdd> (C(D') + K) , (5.82) 

D> D' 

namely that the exponent of the scattering amplitude is invariant under transformations 
of the colour factors 

C(D') -»• C(D') + K (5.83) 

for a given set, where K does not depend on which member (D') of the set one considers. 
Another way of saying this is that any contribution to each colour factor which is inde- 
pendent of the ordering of the gluon attachments on each parton line does not contribute 

14 A systematic way to examine presence of subdivergences is presented in section 6. 
15 Note that this is not true for the trivial case in which R is a 1 x 1 matrix. As discussed in section 3, 
such diagrams have C(D) = C(D), and thus R = 1. 
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to the exponent. This is the multiloop generalisation of a property already noted at two 
loops - see eq. (5.11) - which we may illustrate as follows. Considering a line i with two 
gluon attachments, these give a contribution to the colour factor 

Tf Tf = \ {{Tf, Tf } + [Tf, Tf ]) , (5.84) 

where Tf is a colour generator on line i with adjoint index A, and on the right-hand-side 
we have written the pair of colour matrices as a sum of symmetric and antisymmetric 
parts 16 . The symmetric part does not contribute to the exponent, which we view here 
as a consequence of the symmetry of the exponent under transformations of the form of 
eq. (5.83). That is, the symmetric part of the colour factor does not depend on the ordering 
of the gluon attachments. At higher loop orders, a string of three or more colour matrices 
may occur on any given line, and the above noted property of the mixing matrices again 
corresponds to the fact that the symmetric part of the colour factor does not contribute to 
the exponent of the scattering amplitude. 

In fact, there is more to say about the properties of the mixing matrix in terms of the 
cancellation of sub divergences in the exponent. This is the subject of the following section. 

6. Cancellation of subdivergences 

In the previous section we have seen that the exponent of the eikonal scattering amplitude 
is a sum of webs, where each web (5.79) represents a closed set of diagrams related by 
permuting gluon attachments on the eikonal lines. The mixing matrices of each web have 
zero-sum rows (5.81), which one may understand as a consequence of the vanishing of the 
symmetric parts of the exponentiated colour factors. In this section, we explore these ma- 
trices from a different point of view, namely in terms of the cancellation of subdivergences 
in the exponent. 

As explained in section 2, after performing running-coupling renormalization and ab- 
sorbing any collinear divergences into jet factors, all remaining singularities in the exponent 
of an eikonal scattering amplitude must be associated with the renormalization of the multi- 
eikonal vertex. The latter involves, at each order in the coupling, a single counter-term, 
which removes only simple poles, e^ 1 , in dimensional regularisation. This picture should 
be consistent with the fact that the kinematic factors T(D) corresponding to individual 
diagrams D which appear in the exponent do contain, in general, higher-order poles, up to 
0{l/e n ) at re-loops. Thus, there must be a cancellation of all higher-order poles in the ex- 
ponent. In the following section we perform an exploratory investigation of the mechanisms 
by which this cancellation is realised. 

6.1 Renormalization of the multi-eikonal vertex 

The general structure of singularities of a multi-leg eikonal amplitude is determined by the 
fact that it renormalizes multiplicatively. It would be very interesting to understand how 
this structure emerges in term of webs. 

16 One may view this as a generalisation of the "maximally non-abelian" nature of webs discussed in [7], 
due to the fact that the abelian part of the colour factor does not care about the ordering of gluon emissions. 
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As before, we consider a multi-leg eikonal amplitude, corresponding to the soft matrix 
5(e) in (2.2). To simplify the notation we now drop the colour flow index notation as well 
as the first two arguments of this function (the cusp angles and a s ) but keep e, to stress the 
fact that this function is singular. To be concrete we assume that the eikonal lines are all off 
the lightcone (in the case of lightline eikonal lines there are extra collinear singularities that 
need to be absorbed into colour-singlet jet functions before the following considerations can 
be applied). We apply dimensional regularization and assume, in addition, that there is 
an infrared cutoff so we can systematically work with D < 4 (or e > 0). We emphasize 
that in 5(e) running-coupling renormalization was already applied, and thus all remaining 
singularities are related to the renormalization of the multi-leg eikonal vertex. 

Let us now examine the renormalization of the eikonal vertex along the lines of section 3 
in [102]. We introduce a renormalization factor Z for the multi-leg eikonal vertex, which 
is a matrix in colour flow space, such that 

S Icll (n)=S(e)Z~ 1 (e, f i), (6.1) 

where Z~ x absorbs all the singularities of 5(e), making the renormalized eikonal amplitude 
5 re n(/w) finite, and at the same time ^-dependent. The ordering of the factors on the right- 
hand side of eq. (6.1) is important because of their non-commuting nature. Since Z^ 1 
involves only a single counter-term at each order in the coupling, this equation constrains 
the singularity structure 5(e) may have, thus effectively constraining the corresponding 
webs. 

Let us see how these constraints arise. Both the eikonal amplitude 5(e) and the Z 
factor exponentiate: 

5(e) = exp | £ a> W (e) j Z' 1 (e, /x) = exp j £ <C (n) (e, M ) j (6-2) 

Here we have adopted the notation of [102], in which w^ n \e) represents the sum of all 
diagrams occurring in the exponent at C(a™) 17 . Also, Q n > is a counter-term, consisting of 
pure single pole terms (in the minimal subtraction scheme) which removes, at each order, 
any remaining 0(e~ l ) singularity in w^ n \ As argued in [102], owing to the non-trivial 
matrix structure of the eikonal amplitude and the Z factor, the exponent of 5 rcn (/u) in 
(6.1) is not simply the sum of the exponents of 5 and Z~ x . Rather, defining 

«, = 2>>w c = E a "^ (n) ' ( 6 - 3 ) 

n n 

this sum is corrected by an infinite series of commutators according to the Baker-Campbell- 
Hausdorff formula: 

5 rcn (^) = exp{w(e)} exp{C(e,^)} 

= exp + C + \ [w, C] + ^{[w, [w, (}} - [C, [w, C]]) - ^ [C, [w, [w, C]]] + ■ ■ ■ 

(6.4) 



7 In [102], w (n) is defined so as to implicitly include the factor a™. Here we make this explicit. 
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where for better clarity we did not write the arguments of w and ( in the second line. 

Let us first note that if the [w, C] commutators vanish, as happens for example in the 
two-eikonal line case, the counter-terms £ at any given order in the coupling must simply 
be the negative of the 0(1 /e) term in w. This clearly implies that w is strictly free of 
subdivergences: its singularities are just simple poles. This, of course, goes hand in hand 
with the absence of subdivergences in webs for the case of eikonal amplitudes with two 
lines; we return to this point in the following section. 

In the case of many eikonal lines, things are more complicated. Diagrams which 
contribute to the exponent will indeed contain subdivergences. Yet, the counter-terms 
at any order n, which must render the exponent of 5 ren (^) in (6.4) finite, are strictly 0(e _1 ). 
How is this possible? 

We identify two distinct mechanisms which combine to facilitate the necessary cancel- 
lation of singularities: 

1. The mixing matrices Rdd' 111 (5.79) generate particular linear combinations of kine- 
matic functions in which certain subdivergences cancel. This cancellation occurs 
within sets of graphs W- n ^ (webs of order n, where the subscript i refers to a specific 
web) which build up 18 the contribution to the exponent w(e) of the unrenormalized 
eikonal amplitude S(e) at any given order n: 



w(e) = 



EE^ (n) 



E E E n») <S> c(d') . (6.5) 



D,D' 



2. The commutator terms in eq. (6.4) involving lower-order webs and lower-order counter- 
terms generate multiple epsilon poles. 

We emphasize that despite considerable cancellation through the mixing matrices, 
which we shall elucidate in the following subsections, webs are not entirely free of subdi- 
vergences. Indeed, as we shall see, starting at three loops, at least some multiple poles 
do indeed survive in the exponent of the unrenormalized eikonal amplitude (6.5). This 
is where the second mechanism becomes essential. To appreciate its role, let us examine 
eq. (6.4) at the first few orders in a s . First one may define 



S'ren(^) = exp 



E a > 



n w {n) 
s "^ren 



(6.6) 



where t^ren, as the notation suggests, is the renormalised version of collecting all 
diagrams at 0(a"). Each term u>ren is thus finite as e — > 0. Equating the exponent of 
eq. (6.6) with that of eq. (6.4) at successive perturbative orders, one finds 



w. 



(i) 



(D + C (i). 



w. 



(2) = ^) +C ( 2 ) + ^ u; (l), C (l) 



W 



(3) 



W 



(3) ■ A3) 



^(p),c (2) ] + p ) ,c (1) ]) + ^ 



w 



(1) _ C (l) 



(6.7) 
(6.8) 

(6.9) 



18 It looks natural to assume that there are no further cancellations between different webs W\ n ^ and 
W^ n \ although we have not constructed a proof this never occurs. 
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One now sees explicitly that at first order, the counter-term is simply the negative of the 
single pole in w/ 1 ). Evaluating eq. (6.7) at C(e _1 ) (using the fact that wicl is manifestly 
finite as e — > 0), one finds 

{W = - w W\ e -i. (6.10) 

Things become less trivial already at two-loop order. Individual two-loop diagrams do 
indeed have subdivergences. A priori these could either cancel out in webs owing to the 
mixing matrices, or survive in webs, generating 0(e~ 2 ) contributions to in eq. (6.8). 
In fact, the latter possibility can be excluded since C^ 2 ** = 0{e~ l ) and the commutator 
term in (6.8) cannot contribute at 0(e~ 2 ): such a contribution could only come from the 
0(e~ 1 ) in both and w^, but these are proportional to each other owing to (6.10), 
and therefore commute. There will, however, be a contribution from this commutator at 
C(e _1 ), coming from the O(e ) part of [102]. Consequently, the two-loop counter term 
£( 2 ) is not the negative of the single pole of w^ 2 \ but rather, 

C^^-^l^-^^l^o),^] • (6-11) 

We saw that to conform with renormalization, must only have a single pole, despite 
the fact that some of the contributing diagrams have double poles owing to subdivergences. 
We therefore conclude that at two loops subdivergences must cancel entirely due to the 
mixing matrices Rdd 1 - An example is provided by the diagrams (2a) and (2b) of figure 16, 
which form a two-loop web Wn^u- The contribution of these diagrams to the exponent 
of the unrenormalised soft function is 

where we have expressed the results of section 5.2 in terms of the appropriate mixing 
matrix. Multiplying the vector of kinematic factors by the mixing matrix, this in turn can 
be rewritten as 

*™ - \ tet^r - 1 ™ - ™ ™ - ^ • 

(6.13) 

Both the diagrams (2a) and (2b) contain ultraviolet subdivergences, as can be easily seen 
from the fact that one may shrink the innermost gluon to the origin independently of the 
outermost gluon. Thus, J-(2a) ~ F(2b) ~ e -2 . In fact, an explicit calculation shows that 
the e -2 parts of F(2a) and J- (2b) are equal [58,74]. One then finds immediately from the 
combination of kinematic factors appearing in eq. (6.13) that subdivergences cancel due to 
the mixing matrix. 

It is important to remember that we derived the mixing matrices by considering how 
exponentiated colour factors are related to the conventional colour factors of sets of graphs. 
We now see that this same mixing matrix is responsible for conforming with renormalization 
through the cancellation of subdivergences, implying a relationship between the colour 
structure of graphs, and their kinematic information. The two-loop example we analysed 
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is a special case in that full cancellation of subdivergences is achieved by the mixing matrix. 
It is clear that this cancellation occurs within each closed set of diagrams. 

At higher orders, the mixing matrix will only partially cancel subdivergences in general, 
with remaining higher order e poles being mopped up by the commutator terms of eq. (6.4). 
Because these are built out of lower order counter-term and webs, multiple pole terms in 
iy( n )(e) remain highly constrained: according to eq. (6.4) they can be fully determined by 
lower order webs. 

We note in particular, that the cancellation of the leading subdivergence in webs, at 
any given order n, is complete: w^ n \e) never contains an 0{e~ n ) singularity. We have 
already seen above that at two loop order subdivergences cancel within w/ 2 ) itself, so that 
u/ 2 ) is manifestly 0(e~ l ). It follows that at 0(a 3 ) (i.e. in eq. (6.9)), the only potential 
source of e -3 poles is due to the nested commutator term 



w 



(i) 



c (1) , 



w 



(1) £(1) 



This vanishes at 0(e -3 ) however, from eq. (6.10). It follows that must have no e~ 3 
singularities, and thus that the cancellation of the leading subdivergence relies solely on 
the mixing matrices Rdd> at this order. It is straightforward to iterate this argument to 
higher orders. At 0(a") the leading subdivergence is C(e _n ). Due to the cancellation of 
the leading subdivergence in for any j < n — 1, the only potential source of e~ n poles 
in the commutator terms is via maximally nested commutators of the form 



w 



(i) 



w 



(1) 



w 



(1) >C (1) 



(plus similar terms in which any of the factors is replaced by C^)- However, these all 
vanish from eq. (6.10). We have thus proven that at any order n, 



w 



(n) 



|o(e-™) = • 



(6.14) 



so that the cancellation of the leading subdivergence at any order is purely due to the 
mixing matrix, rather than the lower-order counter-terms. 

Three loop order is the first order where multiple poles survive in webs. To see this 
consider the cancellation of the 0(e -2 ) contribution on the r.h.s. of eq. (6.9). We obtain, 
after some algebra, 



^ (3) le>( £ - 2 ) 



1 

12 



(i) 



w 



(i) 



O(e ))C 



(1) 



(6.15) 



which need not vanish, in general. It is interesting to note, however, that this three-loop 
contribution can be deduced directly from the one-loop calculation. This is of course not 
so for the three-loop counter-term, 

1 



C (3) = — ( %< e -i)-±(p%<^,C (2) 



+ 



w 



(2) 



\o(eO),C {1) }) 



1 

12 



(1) 



(6.16) 



which requires both a two-loop and a three-loop calculation to be performed down to the 
O(e ) and C(e _1 ) contributions, respectively. 
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To summarise, one may identify two mechanisms for the cancellation of sub divergences. 
There are cancellations encoded by the mixing matrices Rdd 1 -, as well as ones due to nested 
commutators of counter-terms and lower order webs. For the leading sub divergence, only 
the mixing matrix is relevant, rendering the highest singularity of 0(a™) webs of order e ra_1 . 
The full interplay of these two cancellation mechanisms clearly deserves further study. For 
now, however, we focus on elucidating the role of the mixing matrices and their properties 
in facilitating the necessary cancellations within webs, and postpone full investigation of 
the counter-term structure to future work. 

6.2 Identifying subdivergences 

To make our discussion of subdivergences more precise, we begin by describing a systematic 
approach to characterising ultraviolet divergences in the relevant type of graphs. Let us 
first reconsider the case of two eikonal lines. Then the absence of subdivergences in webs 
is easily proven [56], using the general rule that the superficial degree of divergence D of 
any (sub-) diagram is given by 

D = 4 - E E f( s f + !) - E N ***> ( 6 - 17 ) 

/ i 

where Ef is the number of external lines of type / (bosons or fermions), and Sf = 0, \ for 
bosons or fermions respectively. Also, Ni is the number of vertices of type i, with 

A i = A-di-^2( 8f + l)n i ,f, (6.18) 
/ 

where di is the number of derivatives in vertex i, and riij is the number of fields of type / 
at that vertex. Note that eikonal lines count as fermions, thus have Sf = \. 

The values of Aj for various different vertex types relevant to our analysis are collected 
in table 2. One sees that Aj = for all vertices we will encounter except the coupling of the 



Vertex 


A, 


EEg 





3<7 





4<7 





nE 


4- |n 



Table 2: Vertex factors of eq. (6.18) for various cases. Here E denotes a generic eikonal line, 3g 
and 4g are the three and four gluon vertices respectively, and nE denotes the coupling of L eikonal 
lines by a hard interaction. 

eikonal lines by a hard interaction. Denoting this simply by A, one may simplify eq. (6.17) 
to give 

D = 4-Y J E f ( Sf + l)-A. (6.19) 

/ 
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Now consider the diagram of figure 20, which one may split into subdiagrams by 
drawing a box around parts of the diagram as shown 19 . As can be seen from the figure, a 




subdiagram of the web may consist of the whole web, or subdiagrams in which one or more 
gluons become external. From eq. (6.17), one finds that figure 20(a) has a superficial degree 
of divergence D = 0, indicating the fact that an overall ultraviolet divergence is present. 
The two subgraphs of figures 20(b) and (c) have D = —2 and D = — 1 respectively i.e. 
no subdivergences are present. The generalisation of this analysis to higher-order webs is 
straightforward. The trivial subgraph consisting of the whole web always has D = 0, as 
one expects from the overall ultraviolet divergence. The fact that no subdivergences are 
present follows directly from two-eikonal line irreducibility, the defining characteristic of 
webs: in any given subgraph, one or more gluons will become external such that D < 0. 
This is not the case for reducible graphs, as one may consider either of the distinct pieces 
as a subgraph by itself, which has D = 0. 

Now let us return the L parton case. For a given subdiagram with L eikonal lines and 
n g external gluon lines, we get 

Y,E f (s f + l) = h + n g 
f 1 

while for the hard interaction vertex we have A = 4 — |L. Thus (6.19) yields: 

D = -n g < 0. (6.20) 

It is then clear that all diagrams in the exponent have D = when the whole diagram 
is considered, reflecting the fact that all diagrams have an overall ultraviolet divergence. 
Furthermore, one may again classify diagrams as reducible or irreducible, where the former 
have colour factors which may be decomposed as C(G) = C(Gi)C(G2) for some G\, Gi- 
Irreducible diagrams have no subdivergences by a similar argument to the two-eikonal line 
case. That is, any box which selects a subgraph of an irreducible diagram will make at 
least one gluon external, such that the degree of divergence becomes negative by eq. (6.20). 

19 We consider here only boxes which include the multi-eikonal vertex. Subdivergences in boxes which do 
not include this vertex are dealt with by renormalisation of the strong coupling constant. 
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6.3 The mixing matrix as a projection operator 

Reducible diagrams will indeed have subdivergences. A full classification of their structure 
requires detailed investigation of their kinematic parts, and is beyond the scope of this 
paper. However, as we have seen already in section 6.1, cancellation of subdivergences 
involves an interplay between commutators of counter-terms and lower order webs, and 
the mixing matrices Rdd 1 - It is thus useful to examine some of the properties of the 
latter, which appear to be completely general 20 . We then interpret them in the context of 
cancellation of subdivergences, and demonstrate this by examining a couple of examples. 
Finally we complete the picture by discussing the related combinatorics of colour factors 
subject to (5.81). 

First, we observe that the mixing matrices R are idempotent, namely, they satisfy 

R 2 = R, (6.21) 

strongly suggesting that they are projection operators. This is indeed the case, as we 
explain below. The idempotence property of R immediately implies that: 

• R is diagonalizable. 

• R has only two distinct eigenvalues: 1 and 0, with degeneracies r and (d — r) respec- 
tively, where d and r denote the dimension and rank of R. 

One may then construct the diagonalising matrix Y whose rows are the left-eigenvectors 
v T of R (or, equivalently, v are the right-eigenvectors of i? T ). That is: 

v T R = v T \ n - YRY- 1 = diag(Ai,A 2 ,...,A d ), (6.22) 

where X n = 1 for n = 1, 2, . . . , r, and A n = for n = (r + 1), . . . , d. 

Let us consider now the vector of kinematic functions T(D) corresponding to the set 
of diagrams {D} in a given web. We may denote this vector by J 7 , and express it via the 
linear combinations formed by the diagonalising matrix, J 7 ^ = Y, in such a way 

that the web in (5.79) can be written as: 

w {m,n2,...,n L ) = F T C = F T RC 

= (j^Y' 1 ^ YRY' 1 (Yc) 

= (j- T y- 1 )diag(A 1 ,A 2 ,...,A d )(yc) (6-23) 

H=l 

From the third line it is clear that owing to the presence of zero eigenvalues, certain entries 
in the vector J-~ T Y^ 1 , containing particular linear combinations of the kinematic functions 

20 These properties have been observed in all examples considered throughout the paper (several sets of 
two, three and four-loop of diagrams). We will not prove that they apply for all mixing matrices at any 
order, although we expect they do. 
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J-(D), do not enter the expression for W, i.e. do not contribute to the exponent. As 
explicitly stated in the last line, only the first r components of T Y _1 , corresponding to 
unit eigenvalues, enter the exponent. 

It should be noted that owing to the degeneracy of the eigenvalues there is some freedom 
in choosing the particular basis of eigenvectors: any linear combination of eigenvectors with 
eigenvalue 1 (0) is also an eigenvector with eigenvalue 1 (0). 

Before proceeding it is worthwhile clarifying the precise sense in which R is a projection 
operator. We wish to explain the physical interpretation of the eigenvalues 1 and 0, and 
the fact that these are the only ones that appear. To this end it is useful to rewrite the 
expression in the last line in (6.23) as follows 

d 

1 — 9(D > r) 

(6.24) 

D=r+1 

where we used the fact that upon summing over all D = 1,2, ... ,d one recovers the unit 
matrix ^d=i ^gd ^dh = ^gh- The first term in the last line in (6.24) is easily recognised 
as the full non-exponentiated result for this class of diagrams. It becomes clear then that 
what is being subtracted from it, i.e. the sum over the subspace corresponding to the 
zero eigenvalue of R, are all the contributions that are generated by the exponentiation of 
lower-order graphs. This is analogous to how eq. (5.26) was used in the previous section to 
extract the exponentiated contributions B n at any given order n, by subtracting the result 
of exponentiating the lower orders from A n . The conclusion is then that the zero eigenvalue 
is associated with contributions that are generated by the exponentiation of lower orders, 
while the unit eigenvalue is associated with new contributions to the exponent. This is 
how R projects out the subspace of contributions which appear in the exponent. 

6.4 Cancellation of subdivergences: three-loop examples 

In the previous section we have seen that the mixing matrices R that define webs in multi- 
leg eikonal amplitudes are very special in that they are idempotent and act as projection 
operators on the vector of kinematic functions (or, equivalently, on the vector of colour 
factors) corresponding to the different diagrams in the web. In this section, we make the 
connection between this property and the cancellation of subdivergences. To this end we 
consider again several sets of diagrams, which are sufficiently complex to feature these 
cancellations, and yet simple enough to facilitate explicit calculations. 

First of all, we wish to illustrate how the leading subdivergence i.e. 0(e~ n ) at O(o") 
cancels out within webs, a property we deduced based on the renormalization structure 
of the eikonal amplitudes in section 6.1 (see eq. (6.14) there). More generally, we expect 
that the particular linear combinations of kinematic functions, the J 7 ^ Y^ 1 entries 21 cor- 
responding to the eigenvalue 1, would have a singularity structure that is consistent with 

21 In contrast, no cancellations are required in the (d — r) entries in T T Y^ 1 which correspond to the 
eigenvalue 0, as these do not contribute to the exponent. 



w, 
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(6.4), where all multiple poles are given by commutators between counter-terms and lower 
order webs. 

The first indication that the linear combinations of kinematic functions correspond- 
ing to the unit eigenvalue work to remove subdivergences is that diagrams which are by 
themselves subdivergence-free reside in this subspace. In other words considering a set of 
diagrams containing a particular diagram D$ that has no subdivergences, we note that a 
basis of eigenvectors with eigenvalue 1 may be chosen such that one of the first r columns of 
Y^ 1 would have a single non-zero entry in the row corresponding to the subdivergence-free 
diagram, so as to pick the corresponding kinematic function F{Dq) alone: 



(F T Y- l ) D = J-(A)) 



3D I D < r . 



(6.25) 



For example considering figure 19 with the mixing matrix Rdd' of eq. (5.80), where the 
eigenvalue 1 appears with degeneracy five, we find that two of the five linear combina- 
tions (J- T I" -1 )/) become trivial, and correspond to selecting one of the subdivergence-free 
staircase diagrams, namely 2], [3, 1], [4, 3], [2, 4]]) or 2], [2, 3], [3, 4], [4, 1]]), re- 

spectively. We conclude that in general some of the entries in (J-~ T Y~ 1 )e) with D < r can 
be immediately recognised to be subdivergence-free, since they correspond to picking in- 
dividual subdivergence-free diagrams, whereas in other entries genuine cancellations must 
take place. 

To observe these cancellations let us reconsider first the example of the web 
corresponding to the diagrams of figure 13. These are all reducible, and contain subdi- 
vergences. By considering sequential shrinking of the gluons towards the hard vertex, or 
alternatively by applying eq. (6.20) to subdiagrams, one finds that each and every one of 
these four diagrams has two subdivergences in addition to the overall divergence, making 
them C(e~ 3 ) in total. According to section 6.1 this maximally-divergent contribution must 
cancel out entirely within the web. 

The mixing matrix in this case is given by eq. (5.77), and its eigenvalues are 



{A}(i,2,2,i) = {1,0,0,0}, 



(6.26) 



so here the dimension and the rank are, respectively, d = 4 and r = 1. Determining the 
eigenvectors we get the diagonalising matrix, which forms the following representation of 
the vector of kinematic factors J-'(D) for this set of diagrams: 

F(D) = Y J HG)Y G l E Y ED 

G,E 

( T(3a) - 2T(3b) - 2T(3c) + T(3d) \ 
1 -T(3a) + 2T(3b) + 2T(3c) + hT(3d) 
~ 6 T(3a) - 2F{3b) + AjF(3c) + F{3d) 
v F{3a) + 4F(36) - 2F{3c) + F{3d) J 



( 


1 


-1 






-1 





1 




2 





1 


\ 


2 


1 


Oy 



(6.27) 



Here we have explicitly written the matrix of left-eigenvectors Y in the second line, and 
absorbed the inverse matrix Y^ 1 into the vector of kinematic parts. Note that the first row 
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of the eigenvector matrix corresponds to the eigenvalue A = 1; all other rows have A = 0. 
As we have seen in eq. (6.23), this tells us that, after acting with Rdd' on the right-hand 
side, only the first kinematic combination in eq. (6.27) contributes to the exponent of the 
scattering amplitude. The other combinations do not contribute to ^1^,2,1) > because they 
are associated with zero eigenvalues of the mixing matrix. Note that we have already seen 
this in eq. (5.16), where the very same combination of kinematic functions corresponding 
to the uppermost entry of eq. (6.27) appears. We have now identified it as the unique 
combination corresponding to a non-vanishing eigenvalue of the mixing matrix. 

As discussed above, all the diagrams in figure 13 are expected to have a leading 0(e -3 ) 
singularity. However, the coefficient of this leading pole should be different in the four 
cases. In diagram (3a), one may shrink the innermost gluon to a point, merging it with the 
cusp at the hard vertex. There are then two ways in which one can shrink the remaining 
gluons, given that they connect different lines. Likewise, in diagram (3d) one may start by 
shrinking either of the two innermost gluons before shrinking the right-hand gluon. Thus, 
there are again two ways of forming the leading e pole. Things are different for diagrams 
(3b) and (3c), where the order of shrinking the gluons is prescribed by the structure of 
the diagram. In (3b), for instance, one must shrink the upper, right-hand and lower 
gluons in that order. Thus, (3b) and (3c) have only one way of forming sub divergences. 
These combinatoric factors are precisely such that the leading epsilon pole cancels in the 
uppermost entry of eq. (6.27). 

The above argument indicates that the leading subdivergence indeed cancels due to 
the properties of the mixing matrix, provided the sub divergences are related by the simple 
combinatoric factors derived above by considering the sequential shrinking of gluons. In 
fact, this reasoning is expected to work for the leading pole, but not beyond this. To see 
this in more detail, let us consider the web VF(2,3,i) corresponding to the six diagrams in 
figure 14. This is a more complicated set of diagrams than those of figure 13, in that 
not all diagrams have a leading subdivergence. Diagram (3B) has no sub divergences, 
just an overall C(e _1 ) singularity; diagrams (3A) and (3C) have one subdivergence, i.e. 
a leading 0(e~ 2 ) singularity; and the last three diagrams, (3D) through (3F), have two 
subdivergences i.e. a leading 0(e -3 ) singularity. 

The eigenvalues of the mixing matrix appearing in eq. (5.78) are 



{A} (2 , 3 ,i) = {1,1,1,0,0,0}. (6-28) 



We again see that the only eigenvalues are 1 and 0, although here both occur with a 
nontrivial degeneracy. In this case the dimension and the rank of R are, respectively, d = 6 
and r = 3. Constructing the matrix of left eigenvectors, the analogue of eq. (6.27) in this 
case is: 
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^ T 4 



/ 4F(3A) + F(3B) - 2J r (3C) + 7"(3D) - 2T(3E) + F(3F) 
-2T(3A) - 2T(W) - 2F{3C) - 2T(3D) + 4T(3E) - 2T(3F) 

-3T(3A) - 3^(35) + 3J"(3C) 
-4T(3A) - F(3B) + 2F(3C) - T{3D) + 2T{3E) + 5T(3F) 
2T(3A) + 2T(3B) + 2T(3C) + 2T{3D) + 2F{3E) + 2T{3F) 
2T(3A) - T(3B) - 4F(3C) + 5T(3D) + 2T(3E) - F(3F) 



T 



2 -2 
1 -1 
1 -2 



1 

-i 



-1 1\ 
-1 1 



1 

1 

10 0/ 

(6.29) 



As before, the matrix on the right-hand side is the diagonalizing matrix Y introduced in 
(6.22). Its rows are the left eigenvectors of R, and we have absorbed the inverse matrix Y^ 1 
into the vector of kinematic functions, forming the linear combinations J 7 ^ Y~ l which are 
explicitly displayed here. The ordering is such that the upper three entries of this vector 
correspond to kinematic combinations with eigenvalue A = 1, and the lower three entries 
to A = 0. 

The structure of sub divergences for these diagrams is not as straightforward as in 
the case of figure 13, due principally to the fact that not every diagram has a leading 
divergence. However, we may easily see that the leading e poles cancel. These come only 
from diagrams (3D) through (3F), where there are no crossed gluons. Furthermore, one 
expects the coefficient of the highest e singularity, 0(e~ 3 ), to be the same in each of these 
three diagrams: in each of them the order of shrinking the gluons towards the hard vertex 
is uniquely prescribed by the structure of the diagram. This is also verified by an explicit 
calculation [115]. 

One may then check that the upper three entries of the kinematic vector in eq. (6.29), 
the ones corresponding to unit eigenvalue, have no leading 0(e~ 3 ) singularities, as these 
cancel in the combination J 7 (3D) — 2J 7 (3E) +J T (3F), whereas this is not true for the lower 
three rows. The latter evidently contain 0(e -3 ) singularities, but because they correspond 
to a vanishing eigenvalue of the mixing matrix, they do not contribute to the W( 2 ,3,i) web. 
As expected, the exponent remains free of such singular terms. 

It is instructive to examine what happens with subleading divergences in this three- 
loop web. Writing the e poles of each kinematic factor as 



m=3 



no) = E -rJh 



(-m) 



(6.30) 



m=l 



we may write the higher pole terms (m = 2, 3, where the m = 3 terms must cancel within 
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the web) in the upper three rows Fj of eq. (6.29) as: 
p _ 1 />(-3) _ ot-(-3) , 7r(-3)\ 

+ ^ {4d 2) - 24e 2) + 4f 2) + - 2^fe 2) ) ; (6-31) 

F 2 = i(-2^- 3 ) + 4^- 3 )-2^- 3 )) 

+ ^ (" 2 ^D 2) + - 2^ F 2) - 2^~ 2) - 2^~ 2) ) ; (6.32) 

F 3 = i(-3^ 2) + 3^ c 2) ). (6.33) 

We have already argued that one expects J 7 ^^ = J 7 ^ 3 ^ = J 7 ^ 3 ^, so the C(e~ 3 ) poles do 
indeed vanish in eqs. (6.31) and (6.32). Considering the subleading (e~ 2 ) poles, an explicit 
calculation [115] 22 reveals that the 0(e~ 2 ) poles of diagrams (3A) and (3C) are related by 

4~A 2) = 2^fe 2) - (6-34) 

Thus, we see explicitly, for example, from eq. (6.33) that the subleading subdivergence 
does not cancel due to the mixing matrix alone: an 0(e~ 2 ) singularity survives in the web 
W(2,3,i)j as indeed required by (6.15). This illustrates the interplay between the web struc- 
ture and the commutator of lower order counter-terms and webs in eq. (6.4) for subleading 
subdivergences. A full exposition of this interplay is beyond the scope of this paper. 

6.5 Discussion 

Let us then recapitulate what we have learnt, and present our general conjecture. Consid- 
ering the web in (5.79), involving a linear combination of kinematic functions F(D) and 
colour factors C(D') of any of the d diagrams in the set (which are mutually related by 
permutations of the gluon attachments along the eikonal lines) we have identified a very 
interesting structure, which we believe to be completely general. This structure is dictated 
by the fact that the mixing matrix Rdd' is idempotent, having only two eigenvalues, 1 
and 0, with degeneracy r and (d — r), respectively. The unit eigenvalue corresponds to 
new contributions to the exponent, while the zero eigenvalue corresponds to contributions 
that are discarded because they have been already accounted for by the exponentiation of 
lower-order graphs. 

Writing the web in this basis in eq. (6.23) we deduce that the mixing matrix can be 
expressed as follows 23 , 

r 

Rdd> = Y, Y dh y hd>- (6-35) 

H=l 

Using this representation of Rdd 1 in (5-79), and summing over D, we get r independent lin- 
ear combinations (H = 1, 2, . . . , r) of kinematic functions, Y2d -^(-^)^dh- Similarly, sum- 
ming over D' we get corresponding linear combinations of colour factors, J2d' Yhd'C(D'). 

22 We have explicitly verified this calculation [115] using our own methods, and found agreement. 
23 Note that the sum over H in (6.35) extends only over the range 1 to r. Clearly if it had extended to d 
instead one would have obtained the identity matrix J2h=i ^gh Yhd = Sqd- 
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Based on our results, we conjecture that each individual order-n web (closed set of 
order-n diagrams) has a singularity structure which is consistent with the renormalization 
properties of the eikonal vertex: 

• The leading singularities, 0(e~ n ), must conspire to cancel exactly between diagrams, 
such that all linear combinations of kinematic functions which enter the exponent of 
the unrenormalized eikonal amplitude (YId -^(^)^dh corresponding to eigenvalue 1), 
have no 0{e~ n ) singularity. 

• When renormalizing the eikonal amplitude, all surviving multiple-pole terms 0{e~ 3 ) 
with j < n — 1 in these linear combinations exactly cancel the contributions of 
commutators of lower-order terms as dictated by eq. (6.4). 

The renormalization structure of the multi-leg eikonal vertex is highly constrained, implying 
that the mixing matrices must encode significant cancellations at any order. It remains 
for future work to determine whether the renormalization procedure can be applied on a 
web-by- web basis, and establish that individual webs do indeed conform with eq. (6.4). 

Our discussion here focused on the kinematic dependence, where appropriate cancel- 
lations must occur. The structure of eq. (6.23) implies that Y acts also on the colour 
factors. Here another general property of this matrix becomes important: each of the 
left-eigenvectors v T corresponding to a unit eigenvalue, i.e. each of the first r rows of the 
diagonalizing matrix Y have entries that sum up to zero, namely 

J2y hd = 0, VJ3" I 1 <H<r. (6.36) 

D 

Note that the zero-sum property does not necessarily hold for those eigenvectors that 
correspond to a zero eigenvalue r + 1 < H < d, as can be checked in the previous examples. 

It is straightforward to see that the zero-sum property in (6.36) is related to the zero- 
sum property of the rows in the matrix R itself, eq. (5.81) above. Indeed, expressing Rod 1 
using Y as in (6.35) and summing over D' , for any fixed D, we have 

d d r r d 

E r dd> =EE Y oh y hd> = E rsh ( E Y *n) = > ( 6 - 37 ) 

D'=l D'=1H=1 H=l D'=l 



where we changed the order of summation and used (6.36), recovering eq. (5.81) above. 

Eq. (6.36) has a clear physical interpretation: according to (6.23) only specific linear 
combinations of colour factors (Y C)h enter the exponent (r such combinations in total, 
H = 1, 2, . . . , r). What characterized these linear combinations is that they are free of any 
component which is independent of the ordering of the gluons along the line, as imposed 
by the symmetry in (5.83), namely 

E Yhd C(D) = £ Y HD (C{D) + K) , (6.38) 

D D 
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where K is D-independent. We thus see that the linear combinations (YC)h (for H = 
1,2, ... ,r) encode the generalization of the antisymmetrisation property familiar from the 
two-loop case (see e.g. eq. (5.11)). 

In conclusion we have seen that webs have a very special structure, summarized by 
eq. (6.23). This implies that the mixing matrices, or equivalently, the corresponding diag- 
onalising matrix Y, is simultaneously responsible for 

• forming linear combinations of kinematic functions of different diagrams, (J 7T Y^ 1 )h, 
that conform with the singularity structure implied by the renormalization of the 
multi-eikonal vertex. In particular, these combinations are free of the leading subdi- 
vergence (eq. (6.14)), and all their lower-order multiple poles correspond to commu- 
tators of lower-order webs and counter-terms. 

• forming corresponding linear combinations of the colour factors (YC)h which are an- 
tisymmetric with respect to permutations of gluons along any of the lines, as required 
by (6.36). 

Thus, these two operations on kinematics and on colour space are intimately related. 

So far we have considered the properties of strictly eikonal gluon emissions or, equiv- 
alently, Wilson lines. In the next section, we discuss how the notion of webs can be gen- 
eralised further to include contributions in scattering amplitudes arising from subleading 
powers in gluon momentum. 



7. Next-to-eikonal webs 



In this section, we consider the application of the results of this paper to the extension of 
resummation methods in multiparton scattering to include corrections beyond the eikonal 
approximation in the momentum of the emitted gluons. 

In more detail, the perturbative expansions of inclusive differential cross-sections (at 
the parton level) have the generic form 



dz 

m,n 



ln m (l - z) 

a nm ( — r ) + v nmS(l -z) + b nm ln m (l - z) + . . . 



l-z 



+ 



(7.1) 



where 1 — z is a dimensionless combination related to the energy carried by soft gluons, 
and the ellipsis denotes terms which are suppressed by powers of (1 — z) and thus non- 
singular as z — > 1. With the exception of collinear singularities, the first set of terms on the 
right-hand-side of eq. (7.1) can be obtained from the eikonal approximation, in which the 
four-momentum of each emitted gluon goes to zero. The "+" prescription and the 6(1 — z) 
terms represent the contribution of purely virtual diagrams. Finally, the last set of terms, 
suppressed by one power of (1 — z) compared to the first set, arises from the next-to-eikonal 
approximation in which gluon momenta may occur to first order in the amplitude. When 
the logarithms in eq. (7.1) become large, fixed-order perturbation theory breaks down and 
one must resum the enhanced terms to all orders in a s . 
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In addition to the large body of work on eikonal resummation, there have also been a 
number of studies focusing on next-to-eikonal (NE) effects. Early approaches included 
subleading effects in the collinear evolution kernel entering the renormalisation group 
equations for the resummation, thus effecting a partial resummation of subeikonal cor- 
rections [116-121]. More recently, efforts have been made to systematically construct NE 
resummation formulae, such as the approach of [82,83], the physical evolution kernel ap- 
proach of [34,35,84-86], and the application of a modified evolution equation for parton 
distributions [87] within the context of threshold resummation [88] . So far none of these 
approaches attempts to classify the diagrams which underpin NE resummation to all orders. 

This was undertaken in [89] , which used path integral methods to rewrite the problem 
of soft gluon resummation in terms of a field theory for the soft gauge field. A set of 
effective Feynman rules was obtained for emission of soft photons or gluons up to next-to- 
eikonal order, and it was shown for the case of scattering involving two partons that up to 
NE order, a given hard scattering amplitude Mo dressed by soft gluons has the following 
schematic form: 



M = Mo exp 



^C(D E )J(D E ) + C(£>ne).F(£>ne) 

£>e One 



(1 + M r ). (7.2) 



Here the D E are eikonal webs, and -Dne are next-to-eikonal webs i.e. two-particle irre- 
ducible subdiagrams which contain one NE Feynman rule, with all other emissions eikonal. 
These formally exponentiate, and have exponentiated colour factors as do the eikonal webs. 
The term M r is a remainder term which does not exponentiate, but which has an iterative 
structure to all orders in perturbation theory. It collects contributions from diagrams in 
which an eikonal gluon is emitted from an external line, and lands inside the hard sub- 
amplitude Mo- Such contributions have already arisen in the literature as part of the 
Low-Burnett-Kroll theorem [122, 123], which was generalised by Del Duca [124] to cor- 
rectly include collinear singularities. The presence of next-to-eikonal webs in the case of 
amplitudes involving two coloured particles suggests that once webs have been generalised 
to the multiparton C 8.S C . £i subset of next-to-eikonal contributions exponentiates as well. 

The proof of web exponentiation in the two parton case of [89] used the path integral 
method for eikonal and next-to-eikonal resummation. This involves non-trivial combina- 
torics in the case of non-abelian theories, due to the noncommuting nature of the source 
terms in the field theory obtained for the soft gauge field. Nevertheless, a classification of 
NE webs was possible using a replica trick argument. This suggests that the generalisation 
of the replica trick presented in the present paper can be used to study next-to-eikonal 
contributions to multiparton scattering. That is, the multiparton version of eq. (7.2) takes 
the form: 



M = Mo exp 



(1 + Mr)- 



£ F(De) Rl D , C(D' E ) + F NE (D NE )R%,C(D' NE ) 

£>e,£>e -Dne, .Dne 

(7.3) 

where the first sum in the exponent goes over the eikonal webs of eq. (1.4), and the second 
represents the next-to-eikonal analogue thereof, where each of the diagrams -Dne contains 
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one NE Feynman rule. Here J-(D) denotes the kinematic part of a diagram, Rdd' the 
block-diagonal mixing matrix (where each block corresponds to a closed set of diagrams 
forming a given web), and C(D') the conventional colour factor of diagram D' . Given the 
different set of diagrams, the next to eikonal mixing matrix Rj$), is of course different from 
its eikonal counterpart R^) D i- 

In the remainder of this section, we argue that one indeed expects next-to-eikonal 
diagrams in which all gluons are external to the hard interaction to exponentiate. Our aim 
is not to undertake a fully comprehensive calculation of NE logarithms, nor to classify the 
resulting mixing matrices. Rather, we wish to clarify the structure of NE corrections in 
multiparton scattering, in line with the comments made for the two parton case in [89] . 

Let us begin by briefly reviewing the path integral formalism we use. The starting point 
is to consider a scattering amplitude for the production of L final state hard particles, each 
of which may emit further soft radiation. One may then separate the gauge field into hard 
and soft modes. For brevity, and to present the argument in its simplest form, we consider 
explicitly scalar emitting particles. The case of fermionic emitting particles is similar, as 
explained in [89] . In abelian gauge theory, the scattering amplitude has the factorised form 

M( Pl ,..., Pn )= I VA^H(x u ...,x L )]J{p k \(S-is)- 1 \x k }(pl + rn 2 )e iS ^. (7.4) 
J k=i 

Here As is the soft gauge field with action S[.A S ], and H(x±, . . .x£) the hard interaction 
producing the emitting scalar particles at positions X{. The factors (pk\(S — ie)~ 1 \x k ) 
represent propagators for a scalar particle in soft background gauge field - with S denoting 
the quadratic operator for the scalar field in the Lagrangian - sandwiched between states 
of given initial position (the points x k at which the particles are created by the hard 
interaction) and given final momentum p k . Both x k and p k are 4- vectors. The explicit 
factors of {p\ + m 2 ) in eq. (7.4) truncate the free propagators associated with the external 
legs, and there is an implicit integration over the positions x k . The propagator factors in 
eq. (7.4) can be represented as first-quantised path integrals [125, 126]. One finds 

(Pfc + m 2 ){ Pk \{S - \eY l \x k ) = e-™- x "f k (oo), (7.5) 



where 



/fc(oo) = f T>y k exp if dt \\y\ + (p k + y k ) ■ A s (x k + p k t + y k (t)) 
Jy k (0)=0 I JO \ z 

+ \d-A s (xi+p f t + y k {t)^j 



(7.6) 



Here y k is the fluctuation about the classical path associated with the scalar particle k, and 
the path integral is over all such fluctuations subject to the boundary conditions of given 
initial position x k and final momentum p k . Substituting this result into eq. (7.4) yields 



M( Pl ,..., PL ) = JvA» Jv yi ...J Vy L H( Xl , . . . ,x L )e iS ^e-^-P 1+ --- +XL - p ^ 



IJ exp i dt Qy^ + ( Pk + y k ) ■ A s + l -d ■ A s ^j 



x 

k 



(7.7) 
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We see that the scattering amplitude has taken the form of a generating functional for 
a quantum field theory for the soft gauge field Ag. The exponent in the second line 
contains terms linear in Ag, which act as sources for A%. These sources are located on the 
external lines, so that the path integral over Ag generates all subdiagrams that connect the 
sources on the external lines, which can be connected or disconnected (they are subdiagrams 
because they do not contain the emitting particles). Thus formulated, the exponentiation 
of soft photon corrections (in terms of connected subdiagrams) is precisely equivalent to 
the well-known exponentiation of connected diagrams in quantum field theory. 

The advantage of the above approach stems from its clear physical interpretation in 
terms of the worldline trajectories of the emitting particles. To see which diagrams expo- 
nentiate, one must calculate the soft gauge field Feynman rules that result after carrying 
out the path integrations over y k in eq. (7.7). This can be done by systematically expanding 
about the classical straight-line trajectory of eq. (3.1). Note that the classical trajectory 
corresponds to the eikonal approximation in which the emitting particles do not recoil. We 
introduce a scaling variable A such that 0(A°) and 0(A _1 ) constitute the eikonal and the 
next-to-eikonal approximations, respectively. 

Up to NE accuracy we can then write eq. (7.7) as 



Mb 1 ...b„ipi,---,PL) = J d d X!...J d d x L H ai ,„ an (x 1 ,...,x L ) 



J k 



(7.8) 



where 24 

f k (A»,[3 k )=Ve W \ig s 



dt 



f3 k ■ A s (f3 k t) + • M(3 k t) + ±p kli nAP(l3 k t) 



~9s 



A 



oo roo 



JO 



dt / dt' 



5(t-t')g^A»(f3 k t)A»(l3 k t') 



+ 9{t - fi kil g av [0M£(&t)] [A v s {(3 k t')\ 

+ imin(t, /V Pku 9 ar [d a A»(P k t)] [d T A»((3 k t')] 



(7.9) 



Here we have explicitly written the integrals over the positions x k , k = 1 . . . L. As explained 
in [89], these can be carried out by expanding the hard interaction and e~ 1Xk ' Pk factors 
about x k = 0, and lead to next-to-eikonal contributions which depend only on the quantum 
numbers of a single parton leg. That is, such contributions do not involve non-trivial colour 
flows between the external legs, thus may be ignored in the following discussion. 

One recognises the first term in the exponent as the Wilson line exponent of eq. (3.2). 
The other terms, by analogy with the eikonal term, act as source terms for the soft gauge 
field. They generate next-to-eikonal Feynman rules that couple the soft gauge field to the 



4 For the explicit derivation of this result, we refer the reader to appendix B of [89]. 
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outgoing parton lines. Carrying out the path integral in eq. (7.8) generates all soft gluon 
subdiagrams which connect the external lines, where each diagram contains at most one 
NE Feynman rule. We do not write down the NE Feyman rules here, but refer the reader 
to [89] for more details. Here we merely note that the NE Feynman rules contain both 
one and two gluon vertices. Furthermore, the precise form of the vertices depends upon 
the spin of the external lines (as is expected for corrections to the eikonal approximation, 
which is itself insensitive to spin effects). The case of scalar particles is shown above. 



Figure 21: Examples of diagrams entering the exponent at next-to-eikonal order, in which • 
denotes a NE Feynman rule. 



The power of the path integral approach used to derive eq. (7.8) is that the argument 
for exponentiation of the scattering amplitude in terms of webs is the same at both eikonal 
and NE order. By analogy with equation (3.5), one may define a generating functional 



Z NE = J [VA»\ e iS ^ fW ® f( 2 ) ® f( 3 ) . . . ® fW, 



(7.10) 



where f( k ^ is the generalised Wilson line operator defined in eq. (7.9), and there are L 
parton lines. Constructing again a replicated theory with N identical copies of the soft 
gauge field, one has 

(Z^f = jvAl...A» « (ff )tf ) . . . f W) ® (ff ) . . . ff) ® 

...® (fW.-.f^) , (7.11) 

(k) 

where is the generalised Wilson line operator associated with parton k and replica 
number i. Each parton line carries a product of generalised Wilson lines. Using the replica 
ordering operator 1Z, such a product may be written as 



-(fc)-(fc) 



f (fc) 



N 



1ZV exp < \g s \ ^ I dt 
{ i= i J ° 

r»oo />oo 



oo 



• N poo poo r -| 

" 9s J E I dt I dt ' [2^ - ^ A^ k t)Al((3 k t>) 



+ ^min(M')/f 4V r [d a Ai{l3 k t)] [d T Ai(P k t')] 



(7.12) 
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Note in particular that in the two gluon vertex terms, the gluons emerging from the 
vertex must have the same replica number. Diagrams in the replicated theory are similar 
to that obtained in the eikonal case. That is, their kinematic parts are the same as in 
the unreplicated theory, but the colour factors correspond to the replica ordered diagrams, 
rather than the original diagrams. One may again write 

(Z NE ) N = l + NlnZ NE + 0(N 2 ), (7.13) 

so that it remains true that diagrams which are linear in the replica number contribute 
to the exponent of the scattering amplitude. Finally, one extracts the coefficient of N 1 
as before and derives the exponentiation of the scattering amplitude up to next-to-eikonal 
order. The exponentiated colour factors are again given by the colour factors of graphs in 
the replicated theory, where the modification is understood as arising from the action of 
the 1Z operator. Examples of next-to-eikonal diagrams which enter at two loop order are 
shown in figure 21. The first two examples enter singly in the exponent i.e. do not mix 
with other diagrams. The third example mixes with similar diagrams in which a different 
eikonal coupling is replaced by the next-to-eikonal vertex. 

As already discussed for the case of two eikonal lines, NE webs are not the only 
source of next-to-eikonal corrections. There are also diagrams in which an eikonal gluon is 
emitted from an external line, and lands inside the hard interaction. These are related to 
derivatives of the hard interaction and jet functions with no emission, using the appropriate 
formulation of the Low-Burnett-Kroll theorem [122-124]. Nevertheless, the results of this 
section show that a subset of NE corrections indeed formally exponentiates. 

8. Conclusions 

In this paper we have investigated the generalisation of soft-gluon exponentiation in terms 
of webs from two-parton scattering to the L-parton case. We found that the idea of webs 
survives in multiparton scattering, but that there are important differences with respect 
to the two-eikonal-line case. Our main result is that one finds closed sets of diagrams, 
related by permuting gluon attachments along each of the eikonal lines, which mix with 
each other in the exponent. The relevant mixing matrices have zero-sum rows, reflecting 
the fact that symmetric combinations of colour matrices do not contribute to the exponent. 
This is the appropriate generalisation of the similar property in the two-eikonal line case, 
originally known as the maximally non-abelian nature of webs. We note that in contrast 
to the two-line case, the topology of those diagrams which contribute to the exponent is 
not necessarily irreducible, as has been demonstrated here in numerous examples. 

If many reducible diagrams are present in the exponent, what is the point of classifying 
webs? At the very least, being able to calculate the exponentiated colour factor of any 
diagram allows one to work directly with the exponent. This produces a considerable 
simplification in e.g. calculations of soft anomalous dimension matrices for soft-gluon 
resummation. Furthermore, classifying the structure of webs has allowed us to generalise 
the arguments given in [89] for the formal exponentiation of a subclass of next-to-eikonal 
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corrections, namely those which arise from diagrams in which all soft gluon emissions are 
external to the hard interaction. 

Using the replica trick we have established an algorithm by which the exponentiated 
colour factor of any multi-loop (and multi-leg) diagram can be directly computed. We have 
also derived a corresponding closed-form expression (eq. (4.3)) in terms of colour factors 
of subdiagrams. We have demonstrated the application of the method, considering sev- 
eral two, three and four loop examples, where in each case we have identified the closed 
set of diagrams that mix with each other and determined their mixing matrix. Further- 
more, considering two non-trivial three-loop classes of diagrams we reproduced the results 
obtained by the replica trick by explicitly exponentiating lower-order diagrams, and sub- 
tracting these from the full unexponentiated sum of diagrams. This alternative calculation 
required manipulating products of colour factors as well as products of integrals of one- 
and two-loop diagrams, bringing them into the form of recognisable three-loop expressions. 
These examples made clear the drastic simplification afforded by the replica-trick based 
method. 

Beyond introducing a new calculational method, our generalization of the notion of 
webs to the multi-leg case is a conceptual step forward. The conventional understanding 
of exponentiation through evolution equations and anomalous dimensions has been very 
effective in clarifying the general singularity structure of amplitudes in the dimensional 
regularization parameter, in explaining the way running coupling effects can be incorpo- 
rated to all orders, and in constraining the kinematic and colour structure of the exponent 
in the massless case [60,63]. However, in general, in the massive case and beyond the 
planar limit, not much is known about the kinematic structure or the colour structure of 
the exponent in multiparton scattering beyond the two-loop level, the state-of-the-art of 
explicit calculations [33,43,68-75]. It is clear though that the exponent has a much simpler 
structure than the amplitude as a whole, opening a window to the all-order structure of 
perturbation theory. We have shown that webs provide a complementary way to under- 
stand exponentiation in multiparton scattering, revealing intriguing connections between 
kinematics and colour space. 

Comparing the conventional picture, of exponentiation in terms of the soft anomalous 
dimension, to the picture of webs is very insightful. The former tells us that in dimensional 
regularization, after removing collinear and running-coupling singularities, the eikonal am- 
plitude can be renormalized by introducing a single 0(l/e) counter-term at any order n 
in perturbation theory. In contrast, a straightforward diagrammatic analysis reveals the 
presence of multiple sub divergences, leading to much stronger singularities up to C(l/e n ). 
In the two-eikonal-line case, no difficulty arises owing to the fact that only very special dia- 
grams contribute to the exponent, those which have no subdivergences. These are precisely 
the diagrams that have irreducible colour structure, the ones we call webs. 

In the multi-leg case a major difficulty arises: individual diagrams that contribute to 
the exponent present subdivergences. In this paper we showed that the problem is solved 
upon generalizing the notion of webs in the multi-leg case to be closed sets of diagrams 
rather than individual ones. These sets naturally arise upon computing the exponentiated 
colour factor of any given diagram as a linear combination of the colour factors of other 
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diagrams. This gives rise to mixing as described by eq. (5.79), where any given set includes 
all the diagrams that are mutually related by permuting the gluon attachments to the 
eikonal lines. We argue that the compound object that is formed upon adding up the 
contributions of all the diagrams in the set, is the proper generalization of a web: it admits 
a generalised antisymmetry property with respect to permutations (realised through the 
zero-sum row property of eq. (5.81)) and it is expected to have a singularity structure that 
conforms with the renormalization of the multi-leg eikonal amplitude. 

The mechanism by which multiple-pole terms associated with sub divergences cancel 
out was investigated here in section 6, revealing a very interesting picture. A first level of 
cancellation occurs within each individual web contributing to the unrenormalized ampli- 
tude. This cancellation is facilitated by the fact that the mixing matrices are idempotent, 
and therefore have just two eigenvalues 1 and 0. Thus only particular linear combinations of 
kinematic functions, those corresponding to eigenvalue 1, enter the exponent; these combi- 
nations turn out to be very special: for example, for any web of order n, these combinations 
are free of the leading 0{\/e n ) singularity. A second level of cancellation takes place upon 
renormalizing the eikonal amplitude. In the exponent of the renormalised amplitude, all 
remaining multiple-pole terms, 0{\/e 3 ) with j < n — 1 at 0(a") get cancelled by nested 
commutators of lower order webs and lower order counter-terms, such that a single 0{l/e) 
counter-term will suffice. For webs to match this structure is highly non-trivial. In partic- 
ular, all remaining multiple-pole terms in webs must be given by nested commutators of 
lower-order webs. 

It is important to emphasize that this picture is still largely conjectural: first, the 
properties of multi-leg webs have not been proven, but rather been deduced by considering 
certain sets of diagrams. Second, it is not obvious that renormalization of the multi-leg 
eikonal amplitude can indeed be performed on a web-by- web basis. In any case, further 
study is needed to fully expose the interplay between the singularity structure of individual 
multi-leg webs and the renormalization of the corresponding amplitude. 

We observed that the structure that emerges in webs (see eq. (6.23)) presents an in- 
triguing relation between kinematics and colour space: the same diagonalising matrix Y 
which generates antisymmetric combinations of colour factors (owing to the zero-sum prop- 
erty in the rows corresponding to unit eigenvalue) also acts on the kinematic functions to 
remove sub divergences. This beautiful structure clearly calls for an in-depth mathemat- 
ical analysis, which will hopefully establish that the above results are indeed completely 
general, and shed light on the mechanism of cancellation of sub divergences. 

We have obtained our results using the replica trick, albeit a simplified form of this 
idea, which was inspired by statistical physics methods [101]. Hence, as a concluding 
thought, it is interesting to ponder whether there are any other unexplored 25 uses of this 
idea in high energy physics, given the many other similarities that exist between the fields 
of high energy and statistical physics. 



For some existing applications, see [127-131]. 
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Note added: As the first version of the present paper became public, A. Mitov, G. 
Sterman and I. Sung have completed a study of related topics [102]. While section 2 
in [102] provides an alternative way to derive exponentiation in terms of webs, entirely 
consistent with our results in section 4.2, section 3 in [102] addresses the renormalization 
of the multideg eikonal vertex, elucidating the role of commutators of lower order webs and 
counter-terms, which we overlooked in our original preprint. In this (published) version of 
the paper we confirm some of the results of this section, and furthermore use it to clarify 
the emerging picture of the singularity structure of webs. 
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A. Alternative derivation of ECF using eq. (4.3) 

In eqs. (5.16) and (5.23) we have given the exponentiated colour factors for the three 
loop diagrams of figures 13 and 14 respectively, derived using the replica trick argument 
as explained in section 3. The aim of this appendix is to rederive these results using the 
explicit combinatoric formula for exponentiated colour factors given in eq. (4.3). This serves 
to demonstrate the application of this formula, as well as to illustrate how it encapsulates 
the result of the replica trick. 

Consider first the diagrams of figure 13. The decompositions of figure 13(a) are shown 
in figure 8. There are five of them, and eq. (4.3) explicitly gives 

1 x C(3a) - \ [{C(2a), C(lc)} + {C(2e), C(lb)} + {C (2d) , C (la)}] + \ [C (la)C (lb)C (Ic) 
+perms] 

= C(3a) - - [C(3c) + C(3a) + C(3a) + C(3d) + C(3b) + C(3a)} 

+ - [C(3c) + C(3a) + C(3d) + C(3d) + C(3a) + C(3b)\ 
= * [C(3a) - C(3b) - C(3c) + C(3d)} , (A.l) 

where we have used the labels of figures 15, 16 and 17. The result indeed agrees with the 
result of eq. (5.12). 
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Figure 13(b) has decompositions (36), (26, lc), (2e, 16), (2d, la) and (la, 16, lc), where 
we use a notation (gi, . . . ,g n ) labelling the diagrams gi in each decomposition. Applying 
eq. (4.3) gives 



C(36) - l - [{(7(26), C(lc)} + {(7(2e), (7(16)} + {(7(2d), C(la)}] + ^ [(7(la)C(16)(7(lc) 



+perms] 
1 



L (7(3a) + (7(36) + C(3c) - C*(3d)] , (A.2) 



agreeing with eq. (5.13). The corresponding results for figures 13(c) and 13(d) are 



C(3c) - ^ [{C(2a),C(lc)} + {(7(2e), C(16)} + {C(2c), C(la)}} + ^ [C(la)C(16)C(lc) 
+perms] 

= \ [-C(3a) + C(36) + C(3c) - C(3d)] (A.3) 



and 



C(3d) - 1 [{C(26), C(lc)} + {C(2e), C(16)} + {C(2c), C(la)}] + ^ [C(la)C(16)C(lc) 
+perms] 

= - [C(3a) - C(36) - C(3c) + C(3d)\ , (A.4) 
6 



respectively, again agreeing with the previous results of eqs. (5.14) and (5.15). 

One may also consider the diagrams of figure 14. Figure 14(a) has decompositions 
(SA), (2g, 16), (2a, la) (repeated twice) and (la, la, 16). Applying eq. (4.3), taking care 
with the repeated decomposition gives 



C(3A) - \ [{C(2g),C(lb)} + 2{C(2a), C(la)}] + \ [C(la)C(la)C(16) + perms] 
= X - [3C(3A) - 3C(3C) - 2C(3D) - 2C{3E) + 4C(3F)] , (A.5) 



which agrees with the result of eq. (5.17). Equivalent formulae for the other diagrams in 
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figure 14 are 

C(3B) - \ [{C(2g), C(16)} + {C(2a), C(la)} + {C(26), C(la)}] 

+ J [C(la)C(la)C(16) + perms] 
= ^ [-3C(3A) + 6(7(35) - 3C(3C) + C(3 J D) - 2C(3£) + C(3F)] ; (A.6) 
C(3C) - J [{C(2 5 ), C(l&)} + 2{C(26), C(lo)}] + J [C(lo)C(la)C(16) + perms] 
= i [-3C(3A) + 3C(3C) + 4C(3D) - 2C(3E) - 2C(3F)} ; (A.7) 
C(3D) - \ [{C(2f), C(lb)} + 2{C(2a), C(la)}] + J [C(la)C(la)C(16) + perms] 
= ^ [C(3D) - 2C(3£) + C(3F)] ; (A.8) 
C(3£7) - 1 [{C(2f), C(lb)} + {C(2a), C(la)} + {C(2b), C(la)}] 

+ J [C(la)C(la)C(l&) + perms] 
= \ [-C(3D) + 2C(3£) - C(3F)\ ; (A.9) 
C(3F) - \ [{C(2f),C(lb)} + 2{C(26), C(lo)}] + * [C(la)C(la)C(16) + perms] 
= ^ [C(3£>) - 2C(3£) + C(3F)] ; (A.10) 
all of which agree with the results of eqs. (5.18-5.22). 

B. Example applications of the generalised Gatheral formula 

Here, for completeness, we provide some examples of the application of the generalised 
Gatheral formula, eq. (4.15), which expresses the conventional colour factor of a graph G 
in terms of lower order exponentiated colour factors. First, let us consider diagram (3a) of 
figure 13. The decompositions of this diagram are given in figure 8, and applying eq. (4.15) 
gives 

C(3a) = C(3d) + - C(2a)C(lc) + C{\c)C(2a) + C(2e)C(16) + C(lb)C{2e) + C(2d)C(\a) 

+C{la)C{2d)} + - \c{\a)C(\b)C(lc) + permsl , (B.l) 
J 6 L 

where we have used the labels introduced in figures 15, 16 and 17. In obtaining this result, 
we have summed over all decompositions, and used the fact that the multiplicity factor 
^G\{m H } = 1 in all cases. Upon substituting C(la) = C(la) etc. and using the results of 
eqs. (5.11-5.15), one may combine the colour factors on the right-hand side of eq. (B.l) to 
obtain C(3D), thus showing that eqs. (4.3) and (4.15) are consistent. 

As a second example, let us consider diagram (3A) of figure 14. This is an interesting 
case, as it involves nontrivial examples of the multiplicity factors N G ^ mH y. Applying 
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eq. (4.15) for this graph gives 

C{3A) = C(3A) + X - C(2g)C(lb) + C{lb)C(2g) + 2C(2o)C(la) + 2C(la)C(2o) 

+ * C(la)C(la)C(lb) + perms , (B.2) 

where we have used the diagram labels in figures 15 and 18 on the right-hand side. Here the 
three terms correspond to decompositions containing 1, 2 and 3 diagrams respectively. In 
the first and second terms, we have used the fact that the multiplicity factor NQ\s mH \ = 1 
in all cases. In the third term, we have used the fact that N G ^ m[1 y = 2, and also instated 
the combinatoric factor (3!2!) _1 from the prefactor (nlflmfl!)" 1 in eq. (4.15). One may 
show using previously obtained results for the modified colour factors that the right-hand 
side is indeed equal to C(3A), again demonstrating the consistency of eqs. (4.3) and (4.15). 

C. Additional examples of three and four loop webs 

In this appendix, we give additional examples of three and four loop webs. Together with 
the other examples in this paper, this completes the survey of all three loop webs which 
occur in the case when the external particles are massless (we do not discuss here connected 
diagrams with a single attachment to each eikonal line: in this case there is no mixing, 
each such diagram is a web by itself and its exponentiated colour factor is equal to the 
original one) . We also give a couple of examples of sets of graphs containing an eikonal- line 
self-energy, which is present only in the case of massive external particles. Finally we give 
another interesting four-loop example. We adopt the notation for diagrams introduced in 
section 3, and present the results for the mixing matrices Rdd' of eq. (5.79) together with 
the column vector C(D') on which it acts. 

Let us begin with an example with a three gluon vertex, shown in figure 22. There 





[[1,2],[1,1],[2],[]] [[1,2],[2,2],[1],[]] 
Figure 22: Diagrams contributing to eq. (C.l). 



is a doublet of diagrams with (2,2,1) attachments, where by this notation we denote the 
number of gluon attachments on lines 1, 2 and 3 respectively. The mixing matrix is given 
by 

! / 3-3\ /c[[l,2],[l,l],[2],[]]\ 
6 [-3 3] ^[[1,2], [2, 2], [!],[]] J ' 
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where we also show the vector of colour factors on the right-hand side, so as to define the 
ordering of rows in the matrix. 






[[1],[1,2,2],[2],[]] [[1],[2,2,1],[2],[]] [[1],[2,1,2],[2],[]] 

Figure 23: Diagrams contributing to cq. (C.2). 



Next, we consider the diagrams of figure 23 i.e. a triplet with (1,3,1) attachments: 

/ 3-3()\ /C[[l],[l,2,2],[2],[]]\ 

-3 3 C[[l],[2,2,l],[2],[]] (C.2) 
V-3-36/ VC[[1],[2,1,2],[2],[]]/ 




[[1], [2,1,3,2], [3],[]] 




[[1], [2,1,3,3], [2], []] 





[[1], [2,3,1,3], [2], []] 




[[1], [2,3,2,1], [3], []] 





[[1], [1,2,3,3], [2], []] 




[[1], [2,2,3,1], [3], []] 





[[1],[2,2,1,3],[3],[]] 




[[1], [2,3,3,1], [2], []] 




[[1],[1,2,2,3],[3],[]] [[1], [2,1,2,3], [3], []] [[1 ],[1 ,2,3,2], [3], []] [[1 ], [2,3,1 ,2],[3],[]] 

Figure 24: Diagrams contributing to cq. (C.3). 

In figures 24 and 25 we show a couple of three-loop sets of diagrams where one of the 
gluons is attached to the same line at both ends, an eikonal-line self-energy type graph. 
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As remarked above, these diagrams vanish in the massless (lightlike eikonal line) case but 
not so in the massive case. The mixing matrix is 
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[[1], [2,3,2], [3,1], []] [[1], [2,3,3], [2,1], []] [[1],[2,2,3],[3,1 ],[]] 




[[1],[2,3,3],[1,2],[]] [[1],[2,2,3],[1,3],[]] [[1], [2,3,2], [1 ,3],[]] 

Figure 25: Diagrams contributing to cq. (C.4). 



We have used the notation (m, n,2, ■ ■ ■ , n£) throughout the paper for indicating the 
number of gluon attachments on each external line. We emphasize that this notation does 
not uniquely specify a given set of diagrams, as can be seen from the following example of 
a (1,3,2) web involving a self-energy. The relevant diagrams are shown in figure 25, and 
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their topologies are different to those obtained by rotating the diagrams of figure 14 to 
form (1,3,2) diagrams. The mixing matrix in this case is given by 



( 3 -2 -1 2 1 -3\ /C[[l], [2,3,2], [3,1], 

l-l-l 1 C[[l],[2,3,3],[2,l], 

-2 2 2 -2 C[[l], [2,2,3], [3,1], 

-2 2 2 -2 C[[l], [2,3,3], [1,2], 

l-l-l 1 C[[l], [2,2,3], [1,3], 

y_ 3 i 2-1-2 3/ \C[[1], [2,3,2], [1,3], 



(C.4) 



Now we have a couple of examples where all three gluons connect different lines, first 
with (2,2,2) attachments, shown in figure 26: 





[[1,2],[2,3],[3,1],[]] 



[[1,2], [3,2], [3,1], []] 





[[1,2],[1,3],[3,2],[]] 



[[1,2], [3,1], [3,2], []] 






[[1,2],[3,2],[1,3],[]] 



[[1,2],[1,3],[2,3],[]] [[1,2],[2,3],[1,3],[]] 

Figure 26: Diagrams contributing to eq. (C.5). 
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(C.5) 
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[[1], [2, 3,1], [2], [3]] [[1], [2, 3,1], [3], [2]] [[1], [2,1 ,3], [3], [2]] 




[[1], [1,2, 3], [3], [2]] [[1], [2, 1,3], [2], [3]] [[1 ],[1 ,2,3],[2],[3]] 

Figure 27: Diagrams contributing to cq. (C.6). 



and then with (1,3,1,1) attachments, shown in figure 27: 
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2/ 




U[[l],[l,2,3],[2] > [3]]/ 



Finally, we consider the four- loop example of figure 28. The mixing matrix in this case 




[[1,2],[2,1],[1,2],[]] [[1,2],[1,2],[1,2],[]] [[1,2],[2,1],[2,1],[]] [[1 ,2],[1 ,2],[2,1],[]] 

Figure 28: Diagrams contributing to cq. (C.7). 
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is given by: 



1 

24 



/24 -24 \ 

-24 24 
\ -24 24/ 



/C[[l,2],[2,l],[l,2],[]]\ 
C[[l,2],[l,2],[l,2],[]] 
C[[l,2],[2, 1],[2, !],[]] 

\C[[1,2],[1,2],[2,1],[]]/ 



(C.7) 



Note that here the three-eikonal-line reducible diagram [[1,2], [1,2], [1,2], []] does not appear 
in the exponent, while all three others, in which the gluons are entangled, do. 
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